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The effect of rectification and non-normal accommodation on the 
two-factor theory of alternating current stimulation is studied. Recti- 
fication with normal accommodation increases the required amplitude 
of the current, especially for high frequencies. Rectification with non- 
normal accommodation introduces two upper limits, one for very low 
and one for very high frequencies. 


The two-factor theories of N. Rashevsky (1948) and A. V. Hill 
(1933, 1936) account in a general way for the results of experiments 
with alternating currents. However if the principal time constant 
is determined independently, marked deviations occur at frequencies 
above 1000 cycles per second. 

One of the possible causes for this deviation with which this 
paper will deal is rectification (Cole and Guttman, 1941). It was 
discovered that the squid giant nerve fiber will pass a greater amount, 
of cathodal than anodal current for a given potential difference. The’ 
rectification ratio was found to be about 20:1. That is, one can say 
in this case that the electrical resistance of the nerve to an anodal 
current is 20 times greater than that for a cathodal current. For the 
frog sciatic nerve (Guttman, 1944) the rectification factor is about 
1.03. The presence of a resistance due to the interstitial fluids makes 
it possible for the intrinsic rectification factor to be considerably 
larger. At the same time this resistance makes it possible for the 
nerve fibers to be under the influence of a partially rectified current 
even though the impressed current is of constant amplitude. This 
point will be discussed in a subsequent paper. tea 

We will therefore break up the cycles into the positive half cycle 
(with a maximum amplitude of I,) and a negative half cycle (with 
a minimum amplitude of of, where ao equals the factor of rectifica- 
tion). Using Rashevsky’s notation, we have 
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de/dt = KI, sin wt — k(e — €0) | 
dj /dt = MI, sin wt — m (7 — Jo) | 


de/dt a= IO 5 sin wt — k(e = 0) 
dj/dt = Mal, sin wt — m(j — jo) | 


2.0; 


(1) 
ie" 


with « — 7 = 0 as the condition for excitation. The angular frequency 
@ is equal to 2a times the frequency in cycles per second. 

For convenience we will only give the solution of the equations 
for the e’s since the equations for the j’s will be entirely analogous. 
These solutions are: 


e=e,+ KI, [k sin wot — w cos wt]/(w? + kh?) + Ae*, I20; 
2 
eé= e+ Kal,[k sin wt — w cos wt]/(w? + hk?) + A’e*, I< 0. (2) 
Since att = 0,¢=«, therefore by solving for the constant for each 
half cycle, putting in the boundary value for the next half cycle and 
using the method of induction, we obtain, for J = 0, 2pn/m < t < 
(2p + 1)a/o, 


E— & 
KI, : e2Ptk/wo —_ 1 
| k sin wt — weos wt + of (1 — a) e*7/ (—) | (3) 
ae ke ekt/o — J] 
+ wen} : 


where the term within the brackets is obtained by the use of the for- 
mula for the sum of a geometric series. As there will be no excita- 
tion when I < 0, the solution for I < 0 is omitted. 

Since we are only interested in steady state responses, we let 


t > o, and therefore p — ow, in which case we obtain, for J => 0, 
0<t<az/o, 


E— Eo 
KI, ek(r/u-t) 
+" | ksin ot — o coset + 0(1—a) aie ag be (4) 


The corresponding expression for I < 0 is similar to this except that 
the sine and cosine terms are multiplied by a. 


We shall not restrict ourselves to normal accommodation. For 
convenience we shall define 6 such that 


M/m=6£K/k. (5) 
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We wish to find the current of minimum amplitude, J,=T,, , which 
will just produce excitation in the steady state. To find this value, 
which will be a function of the frequency and the parameters of the 
nerve, we first obtain « — 7 from equation (4) and its analogue and 
set it and its derivative equal to zero. The first of these equations 
involves (79 — &). This may be eliminated since the rheobasic cur- 
rent for a constant cathodal stimulation is equal to 


Tz =e (Jo — &)/K [1—B + B(L— m/k) (Bm/k)™*™]. (6) 
We will therefore have to solve for the required J,, from the follow- 
ing two equations: 


Teil — 6 +31 —m/k) (Bin/k) ee] 
k2 Boy 2 
a or if — ee | sin wt 


Ge ae hee 


| ee (7) 
k pm eralio- 
= —_— t+ (a1— ——__— ——_—. 
2) es — Pi [eos « ( a) @ ow? + k2 ekt/o J 
pm em (7/o-t) | 
7 wo? + m2 ert/eo — | | Y 
0= i = eee | cos ot + o| p oo Alb | sin wt 
ork ott ml risk may 72 

(8) 


pm? eM (1/w-t) he? ek(m/w-t) 
wz +m? e™™/o—1 a? + kk? eftF/o — 1 ; 
The above equations can be solved explicitly only for very large 


or very small values of w, the angular frequency. For very large w , 
the above reduces to 


Lpfte 
1—p+B(1—m/k) (Bm/k ym" ia 
{(2V2— (1 —2?) k—(1—a) 28k] /220+ (1—«) (1—8) /a} 


which approaches a constant. For k a 2000 and m = 50 sec, this 
constant is 0.92/(1 — a) (1 — 6) if f is not too small. 
For very small ow, 


km[1— B+ 6(—m/k) (Bm/ky”—"Y a 
Vom — ph)? + (1 By em 


Tin/Tr — 
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which also approaches a constant, 0.9/(1 — #) for the same values 


of k and m as above. 
Graphs have been plotted for some values of a and 8. By the use 


of successive approximations, the value of t which would satisfy 
equation (8) for a given value of w was calculated. This value of 
t was then substituted into equation (7) and the solution for In/Ir 


was found giving I,,/Izg as a function of w. 


Dy ties 
(<=0, b= 1)>] ttRelne 
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FIGURE 1 


The experimental curve (solid, with points, Figure 1; Katz; 
1939a, b) follows the simple two-factor theory equation (broken 
curve) for low frequency (Hill, Katz and Solandt, 1936), but for 500 
cycles/sec. and above it deviates from this equation. From 5000 to 
10,000 cycles/sec. the experimental curve follows the equation I,/Ip 
= w/2k. It may be pointed out that in the experiments by B. Katz 
an amplitude modulated frequency technique permits one to disre- 
gard, to a large extent, the effect of a propagated impulse on the 
excitatory and inhibitory states. A value of k = 2100 sec? was de- 
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termined by independent methods. 

The curve (broken) for the two-factor theory with normal ac- 
commodation but without rectification (a = 1, 6 —1) as mentioned 
above, follows the experimental curve for low frequency, but for high 
frequency it is approximated by the equation In/Ip = 0.9 w/k. For 
these as well as all other calculations k and m were taken to be 2000 
and 50 sec"? respectively. 

The curve (upper) for the two-factor theory with complete rec- 
tification and normal accommodation (a = 0, 6=1) deviates from 
the experimental curve at all points and at high frequency represents 
the equation I,,/Ip = 1.8 w/k. 

The curve (solid) for complete rectification but incomplete ac- 
commodation (a = 0, 6 = 0.85) deviates from the experimental curve 
at most points, but not by so much that a suitable choice of a and 
might not eliminate this deviation. The curve flattens out at a high 
frequency which was not experimentally reached. However, the curve 
bends markedly at w = 5000 whereas the data continue upward. The 
value of 6 was estimated from data by H. O. Parrack (1940). 

The curve (dotted) for no rectification but incomplete accommo- 
dation (a = 1, 6 = 0.85) is not appreciably different from the case 
for which a = 1 and 6 = 1 except for low frequencies in which case 
it approaches the curve for which ao = 0, $8 = 0.85. The case in 
which a = 1, 6 #1 does permit the threshold value of the alternat- 
ing current to be appreciably less than rheobasic (Householder, 1944) 
as is found in some experiments (Coppée, 1936). 

Another test of the theory, besides that of carrying out experi- 
ments at very high and very low frequencies, is to excite with shocks 
at various intervals after the nerve has been stimulated with a sub- 
threshold alternating current. In experiments of B. Katz at 5000 
cycles/sec, an amplitude of 0.95 I, , cut off and replaced by a direct 
current, gave regular responses at 0.85 Ip but few at Ip/2. Accord- 
ing to the present simple formulation the required direct current, if 
the change was made at the most favorable point of the cycle, would 
be 0.05 Im = 0.4 Ip since the experimental value of Im/Tr is about 8 

or this frequency. 
; fie gees tine 7 is allowed to elapse with no current, then we find 
that if w > 500, 7 > 3/k = 1.5 msec, 6=1,/In, Io < In and t* is the 
time during the cycle at which the alternating current 1s shut off, 


(e—§)/(é — Jo) =1 + O[B(1 —@) /a— Bm cos wt*/wje". (11) 


For high frequency (w a 4000 per second) (e — 7)/(& — jo) —1 
= 6 Bp(1—a)e"/n, a constant. That is, the amount of shock needed 
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to excite will be greater than the rheobase Jr by a constant factor in- 
dependent of frequency. If there is no rectification, the constant will 
be zero. 

The effect of a local response (Hodgkins, 1938; Katz, 1939) has 
not been discussed in this paper. B. Katz (1939a, b) has suggested 
that this can account for the discrepancy at high frequencies. It may 
be noted, however, that since the time for disappearance of. the local 
response is short, the values computed from equation (11) should 
agree with the experimental results if the theory is adequate. 

It may be that the discrepancy at high frequencies can be ac- 
counted for by the theory of J. Th. Van der Werff (1948). However, 
he has given the results for the first cycle only, the result being the 
same as that given by the two-factor theories. What the result would 
be for the steady state is not clear because of the mathematical diffi- 
culties involved in the solution of the problem. 

The introduction of rectification with normal accommodation does 
not improve agreement with experiment. However, rectification with 
non-normal accommodation does in some respects improve the agree- 
ment between theory and experiment over a considerable range of fre- 
quencies, though it introduces limits on the threshold current for 
which there is no evidence currently available. 

This work was aided by a grant from the Dr. Wallace C. and Clara 
A. Abbott Memorial Fund of The University of Chicago. 
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A group of individuals is considered in which each individual has 
tendencies to exhibit one or another of two mutually exclusive behaviors. 
Neurobiophysically this may be described in terms of Landahl’s re- 
ciprocally inhibited parallel reaction chains. The spontaneous excita- 
tions ¢, and «, at the central connections of each chain are a measure of 


the “natural” tendency of the individual toward one or the other of the 
two behaviors. According to equations derived by H. D. Landahl, the 
probability of one or the other behavior is determined by the difference 


2, — &,. A population of individuals is considered in which ¢, — &, is 


distributed in some continuous way, and therefore in which the prob- 
ability of a given behavior is distributed continuously between 0 and 1. 
The effect of other individuals exhibiting a given behavior is to increase 
the corresponding « of the individual. Thus behavior of others affects 
the probability for a given behavior of each individual. It is shown 
that the equaticns describing the behavior of the population on the basis 
of this neurobiophysical picture reduce in the first approximation to 
the differential equations which were postulated by the author in his 
previous work on social behavior. 


In a recent paper (Rashevsky, 1949, hereinafter referred to as 
loc. cit.) we outlined a neurobiological theory of social behavior based 
on the consideration of learning certain behavior patterns. We de- 
rived rather complicated expressions for the rate of change of the 
number of individuals from one behavior pattern to another. 


In our recent book (Rashevsky, 1948a, hereinafter referred to 
as H. R.) we postulated some simple equations governing this rate 
of change, and we were able to apply those equations to a number of 
social phenomena. It is the purpose of the present paper to show 
how those formally postulated equations may be approximately de- 
rived from neurobiological considerations. 

To this end we shall consider a somewhat simpler mechanism 
than in loc. cit. Instead of Landahl’s “learning circuit,” we shall con- 
sider Landahl’s simple reciprocal inhibition circuit (Rashevsky, 
1948b, chap. xxiv, hereinafter referred to as M. B.), which is re- 
produced on Figure 1 below. 
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FIGURE 1 


Let «, and «, be the amounts of e factors which correspondingly 
excite the circuit so that reaction R, or reaction R, is produced. 
We shall interpret the connection excitations «, and e. as represent- 
ing the desires of the individual for a reaction R, or R, , respectively. 
The threshold h of the system (M. B. page 409) will be taken equal 
to zero. We shall denote the probability of reaction R, by P, and the 
probability of R. by P.. By evaluating the integrals (7) and (8) on 
page 408 of M. B. we find 


for & > &;3 P, =f e*leres 5 P,=1—P,; (1) 
fore,>@; P.=ferers): P,-1—P,: (2) 


where k has the same meaning as in M. B. 

An individual with «, = «, P,; = P. = 4, will have no pref- 
erence for either reaction R, or R.. If «, > e, then P, > P, and the 
individual shows a preference for R,. As e, — & tends to 0, P, 
tends to 1. As e — e, tends to o, P. tends to 1. Hence each indi- 
vidual may be characterized by his P,, which varies from 0 to 1. 
Let N(P,)dP, be the number of individuals having a P, between P, 
and P, + dP. If N, is the total number of individuals, we have 


f rever, Na (3) 
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When an individual is in contact with another individual who 
exhibits behavior R,, and when the latter influences the former to 
exhibit the same behavior, the situation results in an increased stimu- 
lation of the former individual toward activity R,. That means that 
a quantity ¢, is being added to «, so that now we have e, + $; — & 
instead of «, — e.. Similarly, an influence toward performing reac- 
tion Rk, means an addition of a quantity ¢. to e.. Unless di (resp. ¢2) 
is infinite an individual will always occasionally exhibit reaction R . 
The probability of the individual exhibiting reaction R, now (that is, 
under the influence of another individual who exhibits R,) is, accord- 
ing to equation (2), given by 


Pa =)]— 4 e-k(E1+-G1-E2) 7 (4) 


and tends to 1 as ¢, tends to infinity. It is largely a matter of social 
convention as to what probability for R, will be considered suffi- 
ciently large to say that the individual practically always shows re- 
action R,. Let this arbitrary probability, which is close to 1, be de- 
noted by p. In order to insure behavior R, we must have 


shee pee On ee Vi (5) 
or 
1 
ee, ee) (6) 


Thus ¢, must exceed a threshold value given by the right-hand side 
of (6) in order that this influence should insure behavior RF, . For 
&; — & > 0 we have from (2) 


eo — = — log 2(1— Ph). (7) 

For «. — €, > 0 we have from (1) 
a = log 2 Ps, (8) 

Hence the threshold ¢* is equal to 
x" == log ‘— for e, > &, (9) 


and to 


== log 4 Px(1—P) for & > &. (10) 
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Thus ¢,* is a function $*(P:) of P;. The same holds about 2" » 
the threshold of ¢2. 

When the additional strong enough stimulus ¢, or ¢2 is applied, 
the corresponding reaction is not elicited immediately, but after a 
time 7, which is of the form 


TT 6) aR ce log is ; 
Kk t oi om Kz 
where the constants 7, K,, and K, can be expressed in terms of the 
neurobiophysical constants of the neural pathways involved (M. B., 
chap. xxxiii). Thus it takes 7 seconds for an individual to pass into 
the class of people exhibiting behavior Rf, . 
Very roughly we may set the rate of increase dx/dt of the num- 
ber x of individuals which exhibit behavior FR, as 


(11) 


ae 12 
dios So 


If 7) is very small and ¢ sufficiently large, then we have from (11) 
approximately 


1 K, K, 
r= —— log (1- Ja (18) 
K, 1 K,¢1 
Equations (13) and (12) give 
fk ee aa (14) 
dt ee Kes 


A similar relation holds between the rate dy/dt of increase of 
the number of individuals exhibiting behavior R, and the intensity of 
the influence ¢.. 

The quantity ¢, is itself a function of P, of the influencing indi- 
vidual. An individual influences another one toward a given behav- 
ior in the amount in which he himself exhibits this behavior. An in- 
dividual with P, = P; = 4, who when left to himself exhibits equally 
frequently both behaviors, will not influence others one way or an- 
other. We may, for example, set 


$1 =k (P, — Ps) = «(2P,—1), (15) 


when «x is a purely sociological constant determined by the frequency 
of contacts between individuals, and where P, refers to the influenc- 
ing individual. More generally we may put 


os = Uy, (P;,«). (16) 
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Denote by x(P,) the number of individuals with a P, between 
P, and P, + dP,, exhibiting behavior R,. Then the total #; in a 
given individual which results from the influence of all other indi- 
viduals exhibiting behavior R, , is given by 


b= | tu (P, , x) @(P,) dP, . (17) 


Similar considerations apply to ¢., which measures the influence 
toward behavior R,. Denoting by y(P,) the number of individuals 
with a P, between P, and P, + dP, exhibiting behavior R,, and by 
Uz(P;, x) a function corresponding to u, , we have 


= | Ua (P;, x) y(P,) aP, . (18) 
We also have 
x(P,) + y(P,) =N(P,). (19) 
From (14), (17), and (18) we have 
dx (P;) = 5 a ef 
7 =K | u,(P,, x) 2(P,)dP, —K f tn (Py, x) y(P,)dP;. (20) 


An individual for whom 
bt > if uy (P: » x) N(P,) aP, (21) 


will for practical purposes never exhibit behavior Rf, , because even 
if all the other individuals would exert an influence on him in that 
direction, that influence, as given by the right-hand side of (21), 
would not be sufficient. Hence 


$x (Ps) = f "tts (Px, x) N (P,) dP (22) 


gives us an equation for the determination of the quantity P, such 
that all individuals having a P, above that quantity cannot be in- 
duced to perform reaction R,. Let us denote that quantity by 1— 4,. 
In a similar fashion we find a quantity 4, such that all individuals 
having a P, between 0 and A, cannot be induced to perform reaction 
R,. In the terminology of H. R. the two types of individuals would 
be called actives of type I and actives of type II. The total number 
of actives of type I is given by 


ne Pappas (23) 


Ay 
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The total number of actives of type II is given by 
Ao 
Yo= { N(P,)¢Py. (24) 


All other individuals are, in that terminology, passive ones. 
The number of passive individuals exhibiting behavior FR, is 


v1—-Ay 
— | a(P,)dP,, (25) 
As 
while the number of passives exhibiting behavior R, is given by 
1-Ay 
y= {  y(P,)aPy. (26) 
ed As : 
Since the right-hand side of (20) is a constant independent of 
P, we may integrate both sides with respect to P, , between the limits 
0 and 1, remembering that X, does not change with respect to time 


and therefore that d(X, + X)/dt = dX/dt. We find because of (24) 
and (25) 


ax ol wl 
=k | 1, (Pz, «)@(P,)dP,—K | Us (P2, «)y(P:)aP;. (27) 
Now put 


dy = an appropriate average of wu, (P,) 
in the interval (1 — A, ,1); 


a@ = an appropriate average of u,(P,) 
in the interval (A. ,1— 4); 


(28) 
€o>= an appropriate average of u.(P,) 
in the interval (0, 4.) ; 
é¢ = an appropriate average of wu, (P;) 
in the interval (A,,1—A,); 
and denote | 
a—K a; a=Ka; C= A653 o=K é. (29) 
Using the mean value theorem we may now write equation (27) 
dX 
eke + aX —cY — Yo. (30) 


Equation (2) of chapter iii of H. R. is a special case of (30) if 
we make a = ¢ in the latter. 
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It is, however, important to remark that the coefficients a, a, 
¢, and ¢) in equation (30) are in general functions of time, because 
x(P;,) and y(P,) contain time as a parameter. On the other hand, 
in equation (2) of chapter iii of H. R., the coefficients are assumed 
to be constant. If, however, as a crude approximation we substitute 
for the coefficients @ , @, ¢) and ¢ in equation (30) some appropriate 
time averages, then equation (2) of chapter iii of H. R. is obtained 
as a sort of zeroth approximation to equation (30) and hence to 
equation (27). 

Thus we have in principle expressed the coefficients of influence, 
used previously, the terms of the neurobiophysical parameters, K,, 
K,, k, and of the sociological constant «. 


LITERATURE 


Rashevsky, N. 1948a. Mathematical Theory of Human Relations. Bloomington: 
The Principia Press. ' 

Rashevsky, N. 1948b. Mathematical Biophysics. Revised Edition. Chicago: Uni- 
versity of Chicago Press. 

Rashevsky, N. 1949. “Mathematical Biology of Social Behavior: I.” Bull. Math. 
Biophysics, 11, 105-113. 


BULLETIN OF 
MATHEMATICAL BIOPHYSICS 
VOLUME 11, 1949 


INPUT-OUTPUT PROBLEMS IN SIMPLE 
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The dynamics of simple nerve-ganglion systems is studied. The 
frequency of nerve impulses (action currents) leaving the ganglion is 
expressed as a function of the frequency of impulses entering the gang- 
lion. The problem is varied by considering ganglia having internal asso- 
ciation connections. 


Consider a nerve (Figure 1) containing N axones which “map 
into” a ganglion consisting of an equal number of cell bodies. The 
axones of the ganglionic cells, each receiving a connection from the 
“incoming bundle” form the “outgoing bundle.” 


> 
INCOMING — GANGLION OUTGOING 
BUNDLE N CELL BUNDLE 
N FIBERS BODIES N FIBERS 


FIGURE 1 


At the outset it will be necessary to introduce a parameter re- 
ferred to as the “firing frequency” which has the following meaning. 
Arbitrarily choose a “small” interval of time 6 (say from ¢ — 46 to 
t + 46), then the average number of “action currents” occurring 
in a fibre bundle during the chosen interval of time per unit time 
per axone will be called the firing frequency of that bundle at the 
time ¢. 

Now let 4(¢) (the input) and f(t) (the output) denote the firing 
frequencies of the incoming and outgoing bundles respectively. In 
what follows we will investigate the relation between the input and 
output of a system such as the one described above. 

’ To begin with very little can be said about the relation between 
¢(t) and f(t) unless we have some knowledge about the nature of 


the connections found in the ganglion. 
165 
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For some of the simpler cases (e.g., the case in which the con- 
nections are random and the duration of the synaptic delay is the 
same for each fiber) the relation between ¢(t) and f(t) is rather 
evident. In the example cited above f(t) is equal to the value which 
¢(t) had at a time “one synaptic delay” o in the past, thus: 


f(t) =¢(t—o). (1) 


In general, however, we can expect that the synaptic delays will 
vary from cell to cell and will also depend upon the physiological 
environment of the cell as a function of time. In order to include 
the more general case of a ganglion which is characterized not by 
one synaptic delay but by a distribution of them, we will introduce 
the notation K («) de which denotes the probability that an arbitrarily 
chosen neuron is characterized by a synaptic delay in the range 
o—4dc too + 4dc. 

Under these circumstances we can expect that f(t) will be some 
transformation of ¢(¢) which depends upon the form of K(c). More 
specifically, since ¢(t — c) is the firing frequency of the incoming 
bundle at the time t — oc, the product ¢(t — o)K(c)do is the frac- 
tion of f(t) which is due to the input at the time ¢ — c. It follows 
then that if ¢(¢) = 0 for the time interval t = — » tot = 0 that 


re) = | Keys dee (2) 


By a simple transformation of variables this equation can be 
transformed into a more familiar form. For this purpose let 
o = t— 7 whereupon do = — d-. It follows then that when c= 0, 
¢=7 and wheno—t,7=—0. Equation (1) then becomes 


t 
f(t) = { USariars: (3) 


Of course very little more can be said about the relation between 
input and output than is expressed in equation (2) until more is 
known about the form of the function K(c). In any study of actual 
nerve-ganglion systems the form of K(c) would need to be deter- 
mined empirically. For the purpose of our analysis, however, we 
will assume a form of K(c) in order to determine the “kinds” of 
relations which such a system implies. For this purpose let K(c) =~ 
Ke-*° where K is a constant. Inserting this form of K(c) into equa- 
tion (2) we obtain 


yee eee i eX (7) dr. (4) 


0 
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The integral of equation (3) can be formally evaluated if the 
operation of integration by parts is performed an indefinite number 
. of times. 

Using this device, equation (3) after n integrations can be writ- 
ten as 


f(t) =$(6) rue gee (£5 2S yen AO 0 (t) 
K ke? Km 
—e*| 60) 90) += 9"(0 5 
x? ze (0) (5) 
ie 
—— ee ee n-1 +R; 
fet Pe (0) 


i 
i 
where R,, = (—1)” Ke-*' [ = eXt 6™ (7) dr. 


Equation (5) implies that the output f(t) is no longer simply 
a temporal translation of the input as was the case for a constant 
synaptic delay. On the contrary, when the synaptic delays are dis- 
tributed as indicated above, the output may be a serious distortion 
of the input. 

It follows from equation (5) that for t >> 0 we can write ap- 
proximately 


==] n-1 
wih ‘sate 


o- (t) + Ry. (6 
Beer at) (6) 


: if aL erm 
F(t) =$(t) —-—9'(t) += 9" (t) 


It can be seen from equation (6) that if ¢(¢) is a constant then 
after a sufficiently long time f(t) = ¢(¢). 

Note, however, that if ¢(t) is analytic and o a constant suffi- 
ciently small, we can write by means of Taylor’s theorem 


o(t—o) =) —eg'(t) tog" (t) — + Re. (7) 


By comparing equations (6) and (7), it is evident that if {t) 
is a linear function of time, i.e., if all derivatives after the first are 


zero then 


f(t) =$¢(t—o). (8) 


In other words, if the input is a linear function of time, then 
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after a sufficiently long time the output will be the same linear func- 


— 1 : : : 
tion with a time lag of « = i It is interesting to note at this point 


that o = R is the average synaptic delay of all the cells in the gan- 


glion. 

For polynomials of degree greater than one we might well ex- 
pect more serious distortions. In fact for a parabolic input, the out- 
put is distorted by a constant, and in general for an input repre- 
sented by a polynomial of degree n, the output is distorted by a 
polynomial of degree n — 2. 

There are a great many important factors found in actual nerv- 
ous systems which are completely neglected by this model. Never- 
theless it is interesting to note, if only in passing, that the distor- 
tions of input which were discussed above are suggestive of such 
phenomena as adaptation and facilitation. 

Consider now the case in which the ganglion neurons not only 
receive connections from the incoming bundle but also have an “‘asso- 
ciation mapping.” In other words each neuron in the ganglion re- 
ceives one connection from another ganglion cell and also “‘sends out” 
one connection to another ganglion cell. 

In this case the output f(t) is not due simply to ¢(¢t) and K(c) 
but also is influenced by the self-mapping. 

Before considering the general equation governing the relation 
between the input and output let us examine the implications of a 
system of three neurons connected in the manner depicted by Fig- 
ure 2. 


FIGURE 2 


Let f,, fg and f; be the firing frequencies of the respective neu- 
rons, and assume that the output neuron will fire if either neuron 1 
or neuron 2 or both 1 and 2 fire. The expectancy of firing of neuron 
1 for a picked interval of time 6 is given by 6f, and similarly that of 
neuron 2 is 6f,. Now since the probability that 1 and 2 will fire at 
exactly the same time is zero, it follows that if the refractory period 
of each neuron is also zero then 


Ofs = of, + dfs, (9) 
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whereupon we obtain 
fe=fitf.. (10) 


Equation (10) which tells us that the output frequency is simply 
the sum of the two input frequencies can be considered valid for those 
systems in which the refractory period is vanishingly small. We 
know however that this is usually not the case. Therefore, in order 
to account for those times when neurons 1 and 2 both fire within a 
period of time less than or equal to the refractory period of neuron 
3 (and therefore are counted only as a single impulse) , equation (9) 
must be modified to read 


Of, = Of1 + 6fe — Of rf, Cis) 


where now the symbol 6 refers to the magnitude of the refractory 
period. It follows then that 


fe=fa + fe — Ofsfo- (12) 


The output frequency f; is seen then to be a function not only of 
the input frequencies but also of the refractory period 6. 

Returning now to the nerve-ganglion model with the additional 
assumption that the ganglion has a “self-mapping” we find that each 
outgoing axone is a member of a three neuron system such as the 
one just described. In this model f, corresponds to the input fre- 
quency and f, corresponds to the “self-mapping”’ or output frequency. 

Reasoning as before but this time taking into account the self- 
mapping we obtain 


f(t) = [° K(@)6(t— 0) da Ji J K@it—a do 


— a] [Kyoto de | f K@re- ode 


In formulating equation (13) it was assumed that the self-map- 
ping was complete, i.e., each member of the ganglion possesses a self- 
mapping axone. If we change this requirement to the more general 
case of a ganglion having only a fraction 4 of its neurons in a self- 
mapping configuration we must modify equation (18) slightly. For 
this case we have 


ray = [Ke o(t— ode +4 [ Kofta de 


(13) 


(14) 
#4] [Kept orae| [ Ke@re—ode| 
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Now if we consider as previously that K (c) = Ke*? and ¢(t) =90 
for — 0 <t < 0, and again as before if we put ¢ — «o = 7, then 
equation (14) can be written as 


f(t) =Ke*! [ olnetde taken! ff (syer de 
(15) 


ae al Ke-*t i (7) ex ds | Ke-Kt [toe dz | 


For certain simple forms of the input ¢(¢), equation (15) can 
be readily solved. 

We will here consider only the case of constant input, i.e., where 
¢(t) = ¢) = constant. Furthermore, only the asymptotic behavior 
of f(t) will be studied, that is, f(t) will be examined for large val- 
ues of ¢. 

By differentiating equation (15) and letting 


at 
1= Ker" | f (2) e87 dr 
0 


we obtain 


f(t) = Kg ,e*' 
(16) 
HAL = $6.05 66") I —=4k 6,662 
However from the definition of J it can be seen that 
I'=Kf(t) —KlI. (17) 


By solving equation (15) for J and substituting this form into equa- 
tions (16) and (17) we obtain a system of two equations with two 
unknowns. By simultaneous solution we then obtain f(t) as a func- 
tion of t, 65, K, 6,4, and f(t). 

By considering the structure of the nerve-ganglion model we can 
infer that for a constant input, f(t) must be a monotonically increas- 
ing function. However due to the finite refractory periods of the 
neurons, f(t) has an upper limit of 1/6. From this we can conclude 
that for a sufficiently large t, f(t) is approximately constant. In 
other words f’'(«o) =0. 

Using this information it becomes possible to solve for f(o), 
and obtain 

do 
1—A(1—8 go) 


Note that for 4 = 1, that is, for the case of a nerve-ganglion 


FU 8) p (18) 


ALFONSO SHIMBEL 171 


system with a complete self-mapping, equation (18) reduces to 
f(«) = 1/6 quite regardless of the magnitude of the constant input 
of ¢.. This latter result, however, can be immediately predicted by 
examining the model under consideration when 4= 1. For this case 
every ganglionic cell has a self-mapping axone so that no matter how 
small the constant input $. may be, if we wait sufficiently long all of 
them will eventually be firing at their maximum frequency. As was 
pointed out before, this maximum frequency is limited by the finite 
magnitude of the refractory period and equals the value 1/6. 

Note also that for the other extreme value of 1, namely, 1 = 0, 
the solution is again trivial, that is, f(0«o) = ¢). In other words, if 
the ganglion cells have no self-mapping, then the input goes through 
the ganglion undistorted. 

In general, however, for all 4 such that 0 < 4 < 1 we must refer 
to equation (18). 

The author is indebted to Mr. George Karreman for a critical 
reading of this paper. 

This work was aided by a grant from the the Dr. Wallace C. and 
Clara A. Abbott Memorial Fund of The University of Chicago. 
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A tentative suggestion is made to interpret the phenomenon of 
mitosis as being due to coiling and uncoiling of long chain molecules, 
which are assumed to constitute the spindle. A possible quantitative ap- 
proach is briefly outlined. 


The theory of metabolic diffusion forces has proven to be rather 
useful in the theory of the general phenomenon of cell division, and 
has led to correct predictions of the over-all quantitative features of 
this process (Landahl, 1942; Rashevsky, 1948). While the present 
developments do not by any means prove that the diffusion forces 
are the only, or even the principal, cause of cell division, they dem- 
onstrate definitely that in many instances they are a sufficient cause. 
Moreover the correct prediction of the time dependence of the elon- 
gation and constriction of freely dividing cells shows that even if 
other factors should eventually be found to be the cause of cell divi- 
sion, any new theory would have to lead to the same, or to very simi- 
lar, equations for elongation and constriction. 

The situation is, however, quite unsatisfactory when we consider 
other important special aspects of cell division such as mitosis. There 
is no “typical” mitosis any more than there is a “typical” cell division. 
The process of mitosis in different types of cells differs considerably 
in details. One general feature, however, which is common to all 
cases stands out: Mitosis provides a mechanism for an equal distri- 
bution of chromosomes in the daughter-cells. The arrangement of 
chromosomes in the metaphase in an equatorial plate and the subse- 
quent movement of each set of daughter-chromosomes toward the 
opposite poles of the achromatic spindle is a general feature common 
to all mitoses. The greater complexity of the whole phenomenon, as 
compared to the relative simplicity of the over-all process of cell di- 
vision, makes a simple explanation rather unlikely. The more com- 
plex an explanation we may have to use, the more numerous will be 
such explanations. The only proper way of procedure would be the 
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one used previously, namely, to investigate all, or as many as Pos- 
sible, theoretically conceivable explanations in abstracto, then study 
the mathematical consequences of all those approaches, and finally, 
from comparison of those consequences with observation, decide 
which of the approaches is the most useful. 

Whatever the nature of the chromosomes, if their arrangement 
in an equatorial plate should represent, at least temporarily, a stable 
configuration, there are only two conceivable types of forces that 
would produce such a configuration. Either each chromosome must 
be repelled by the poles with a force which decreases with the dis- 
tance from the center of repulsion, or each chromosome must be at- 
tracted toward each pole by a force which increases with the distance 
from the center of attraction. If the poles are metabolic centers pro- 
ducing some metabolite S the first type of forces could be interpreted 
as diffusion forces. If the chromosomes are attached to each pole by 
contractile fibers whose force of contraction increases with their ex- 
tension, we will have the second type of forces. A combination of 
both forces was suggested in an outline of a mathematical theory of 
mitosis by N. Rashevsky (1938). Some quantitative consequences 
were drawn from the assumption that the daughter-chromosomes are 
pulled apart by contractile fibers, and were compared to the meager 
data available (Rashevsky, 1941). As F. Schrader (1944) pointed 
out, however, the nature of the data used was such as to make the 
apparent agreement meaningless. 

The shortcomings of the approach outlined above do not yet 
necessarily preclude its usefulness if it is developed further and im- 
proved. In the present paper we shall, however, outline a different 
possible approach. 

It is generally accepted that the chromosomes are aggregates of 
molecules which perform important catalytic functions in the life of 
the cell. In a sense each of them may be viewed as a giant molecule 
which consists of complex parts, the genes. Argument about termi- 
nology would be rather fruitless since we may refer to a polymerized 
molecule as one giant molecule or as an aggregate of smaller ones. 

The circumstance that usually the chromosomes are not visible 
in the resting nucleus may be interpreted as follows. The chromo- 
some may be viewed as a long chain molecule whose cross dimensions 
are too small to make it visible. During the process of mitosis the 
action of certain chemicals, or perhaps the change in pH , makes each 
chromosome the center of precipitation of chromatin and thus makes it 
visible. The situation is not unlike making individual ions visible 
in a Wilson cloud-chamber. 
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Though not conclusive there is evidence that each chromosome 
is attached at its spindle attachment to a fiber which forms a part of 
the achromatic spindle. Let us assume that what we observe as a 
chromosome is only the temporarily visible part of a T-shaped mole- 
cule, which we shall denote as a C-molecule, as shown in Figure la. 
The chromosome proper we shall designate as the V¢=part4otethe 
C-molecule. The other part of it we shall call the f-part. Let all 
C-molecules have f-parts of the same or approximately the same 
length. 


FIGURE la FIGURE 1b 


During the life cycle of the cell such a molecule is duplicated by 
the formation of a sort of “mirror image” (Figure 1b) and then split 
into two daughter molecules. It is a problem of molecular physics 
to develop a theory of such self-duplicating molecules. 

As a catalyst each chromosome is responsible for the formation 
of some substance. Naturally the production of each substance will 
depend on the particular “life-phase”’ of the chromosome. Some of 
these substances, S, may be produced mostly during the period fol- 
lowing the formation of the daughter chromosomes. During the time 
just preceding the separation of the daughter chromosomes when the 
C-molecules are in configuration, shown in Figure 1b, the produc- 
tion of S will be at a minimum. 

Due to the thermal agitation the C-molecules will be “knotted” 


AY IM 
W/ 


FIGURE 2a FIGURE 2b 
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at random as discussed by W. Kuhn (1934, 1936). Suppose that at 
the time when all C-molecules have the configuration shown on Fig- 
ure 1b something causes their f-parts to orient themselves parallel 
to each other. A configuration shown in Figure 2a will then result 
which reminds us of the general shape of a metaphase nucleus. The 
reason for such an orientation may be either an external force or a 
sort of “inner molecular field,” which may appear at a proper time. 
The appearance of such an orienting field may be connected with the 
presence or absence of some metabolites. For instance, the presence 
of S may cause the weakening or complete abolition of certain chemi- 
cal bonds which tend to hold the f-parts parallel to each other. In 
that case when the S-substance is at a minimum, these bonds will 
come into play and orient the f-parts. Simple quantitative estimates 
which we shall discuss presently indicate, however, that with plausible 
bond energies there can be no self-orienting effect of the f-parts, if 
the total number of C-molecules is small, that is, of the order of ten 
or a hundred. If, however, we assume the presence of numerous long 
molecules in the nucleus which are not all attached to the chromo- 
somes (f:-molecules), then the decrease in a certain S-substance may 
well cause the orientation of the f,- molecules. This orientation cre- 
ates an “inner molecular field” which also orients the f-parts of the 
C-molecules. Thus the picture shown in Figure 2b results. Now if 
the c-parts of the C- molecules separate, either due to the tension of 
the f-parts or for other reasons, this will result in an increased pro- 
duction of S, which will loosen the bonds between the f,—molecules 
and the f-parts and result in the latter becoming “knotted” again. 
This knotting results in a contraction of the f-chains, and thus pulls 
the c-parts, or daughter chromosomes, away. 


If the f,-chains are also susceptible to the action of S , but to a 
different degree so that they break in the middle after the chromo- 
somes have been pulled apart toward their respective poles, then the 
f,-chains also will break and coil up, forming the daughter nuclei 
with the halves of the f-chains. 


We may make, if we wish, a more specific assumption. Let the 
f- and f,-chains consist of ionized molecules, the ionization depend- 
ing on the pH of the medium. The pH in this case may play the role 
of the metabolite, S. At such values of pH at which the dissociation 
is strongest, the chain will tend to elongate because of the electro- 
static repulsion between neighboring ionized molecules. If the divi- 
sion of chromosomes affects the pH in such a way as to reduce the 
dissociation, the chain will coil up under the influence of the thermal 
agitation which now will prevail over the electrostatic repulsion. 
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A tendency toward parallel orientation may also be produced by 
van der Waalsian attraction forces between the end molecules of 
each chain. 

Let us now see what quantitative conclusions could be drawn 
from the above picture. 

The reader who is familiar with the work on the kinetics of long 
molecules will notice that the following discussion can readily be im- 
proved and refined by the use of methods developed by other authors. 
In particular, the problem of “relaxation time” can be treated rigor- 
ously by the method used by J. G. Kirkwood and R. Fuoss (1941). 
However at the present stage of our development we are interested 
only in very rough approximations. There are practically no quanti- 
tative data on mitosis with which to compare any exact mathematical 
conclusions. Therefore we shall confine ourselves to much simpler 
methods which give us only orders of magnitudes. The outlined treat- 
ment will have to be mathematically refined in future developments. 

Consider a long molecular chain consisting of N links, each of 
length 1. In the absence of any external or internal fields, the mole- 
cule, as shown by W. Kuhn (1934) will have a “knotted”? shape due 
to thermal agitation. The probability W(7r,)dr, that the ends A 
and A’ of the molecular chain are at a distance between 7, and 7, + dry, 
is given by 


1 
bY2aN 


W. (rsdn =( ) e wn dnredr, (1) 
where 
raed | 1 + cos f 
V3 WV 1—cos 


6 denoting the valence angle between two consecutive links. Thus } 


is of the same order of magnitude as /. ay 
The probability that the middle point of the chain is at a dis- 
tance between 7, and 7 + dr. from the line AA’ is 


(2) 


27? 


1 eine 
2 (1s dr, =——e wwredre. (3) 
Ws (2) trey : 


Finally, Kuhn introduces a “third dimension” of the molecule in the 
direction normal to the 7,, 7 plane. The probability in this direc- 


tion is given by 


Ws; (r;) dr; — S e wv drs. (4) 
b 


VNa 
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The average values of 7: , 12 and r,, and the mean squares are given 
by 


ono Sees b 
=e [aNo; fae ne Nes (5 NG (5). 


2 2Va 
— — nse ub? ) 
r? = 38ND?; %2? = 4 Nb’; Ess ae (6) 


If 1 ~ 10-* em, then, in order to have a chain of a length ~ 10° 
cm, which is the size of the achromatic spindle, we must have N ~ 10°. 
For such long chains we may apply with Kuhn the concept of entropy 
to a single chain. Denoting by k the Bolzmann constant, we have for 
the entropy of a chain 


s=s,+s,+8,=—kf[log W.(7,) + log W.(r2) + log W3(7s)]. (7) 


Hence denoting by ¢, , ¢, and ¢; three constants we have 


hire 
= ese 8 
$8, =¢,+ 2k log 7, See (8) 
Diets” 
Sao tag OR Tobe sae (9) 
Akr,? 
3. ¢; = aN (10) 


The total work done in stretching a chain from its average length 


7, to its full length IN is given Py, the corresponding change Af of the 
free energy. Since 


f=u—Ts, (11) 


where u is the total energy and T is the absolute temperature, and 
since we consider now the case in which Au = 0, therefore 


Af =— TAs =— T (As, + As. + Ass). (12) 


With b ~ 10% cm, N ~ 10°, we find that 7, ~ 5 X 10-6, 7%? ~ 3 X 
10-" and, from equation (8), As, ~ — 10™ erg degree. The con- 
tribution of the log r, term is insignificant. As the molecule is 
stretched completely the “cross-dimensions” 7, and 7, change from 
their initial values to about 10-° cm. With the same values as before 
we find that As, ~ As, ~ 10°* erg degree?. Hence in computing 
the work and the forces, we may neglect the terms s, and s; in equa- 
tion (7) and put approximately 
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S=S8,. (13) 
From equation (12) it follows that the force of contraction of a 
stretched chain is given by F = — T (ds/or,) and because of (13) 
and (8) we find 
2 r 
F=—kT ( = — ) 
Tee N: 3 


For sufficiently large extensions (7; > 10-5 cm) the first term 
becomes negligible and the force varies approximately linearly with 
the extension. 

If such a stretched chain is left to itself it will again return, 
due to thermal agitation, to the “knotted” configuration described 
by equations (1), (3), and (4). The rate of return may be approxi- 
mately evaluated in the following way: Consider each link of the 
chain as a sphere of radius a = 10-° cm. While the ends of the chain 
move by a given amount toward each other, the middle link remains 
approximately undisplaced, while the velocity of the other links with 
respect to the medium varies between (1/2) (d7r,/dt) and zero. Hence, 
approximately, we may put the average velocity of each link with re- 
spect to the medium as equal to v = — (1/4) (dr,/dt). Assuming 
Stokes’ law, the average force of resistance on each link equals 627 av 
where 7 is the viscosity of the medium. The total force of resist- 


3 
ance to the N links is then 6 a 7 avN, or — ote aN (dr,/dt). Neg- 


lecting inertial terms, this must be equal to the force F as given by 
equation (14) in which, as we have seen, we may neglect the term 
2/r,. We thus find 


ee (15) 


Integrating and taking r, = IN for t = 0 (completely stretched 
chain) we find 
— IN €- T/31 abn?) t 2 (16) 
For values of a, b and N assumed above, for T ~ 300°, and 
n ~ 1 abs. c.g.s units, the value of the exponent is about 1 sec. Thus 
after one second the chain will have “coiled up” to only 1/e of its 
original length and the “relaxation time” is thus of the order of sec- 
onds. 
For a constant length of the chain b*N? = const. Hence the re- 
laxation time increases with a ~ b. For a ~ 10*% cm, N ~ 10°, and 
n ~ 1 .g.s units, we find for the relaxation time a value of about 
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100 sec. With the same values of a and N,, but with y ~ 10°, it again 
reduces to about 1 second. 

The total work of stretching a ain is, according to equation 
(14), approximately equal to kTI2N?/2b°N ~ kTN/2. The energies 
of molecular interactions which may depend on their relative orien- 
tations are of the order of kT ~ 10-* erg, frequently several times 
that amount. Hence if each link contains a molecule or group of mole- 
cules of a polar nature which tends to orient the corresponding link 
of another chain in a parallel direction, the total decrease of the po- 
tential energy when two chains are stretched and arranged side by 
side parallel to each other may well be sufficient to supply the work 
kTN/2 necessary for stretching each chain. In that case a spon- 
taneous orientation would occur. However, the intermolecular forces 
which tend to produce any orientation are of short range and become 
negligible at distances exceeding 10-7’ em. Hence if we had only 
f-fibers, which would number about 10—100 for the whole nucleus, 
no orienting effect would occur because the average distance between 
such fibers would be about 10-* em, which is much too far away for 
any ordinary intermolecular forces to become effective. As we said 
above, however, we may have a large number of f,-molecules. If 
their number equals or exceeds 10”, the average distance will be of 
the order of 10-° cm, and an orientation may occur. 

The orienting forces may act like a “molecular field,” similar 
to the one introduced by P. Weiss (1907) in the theory of ferromag- 
netism. The more each chain is stretched, the stronger may be the 
orienting force, and that in its turn will tend to stretch the chains 
further. Thus, approximately, we shall have a “stretching force” F, 
acting on each molecule and proportional to 7, (which we may con- 
sider as approximately the same for all molecules). Thus 


F,=— Ari, (17) 


the minus sign indicating that the force is a stretching one. For a 
spontaneous orientation to occur we must have F’,, > F, or 


A— (kT/0?9N) > 0. (18) 


In that case we may roughly estimate the rate of orientation, The 
viscous resistance is again equal to (3/2)a 7 aN (dr,/dt). This must 
be equal to [A — (kT/b?N)]7,. Thus we find now 


o ar; kT 
—n2n aN —=|(A— )n. (19) 


Putting now for t= 0, 7, = 7,, as given by equation (5), and intro- 
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ducing the notation 


2 A — (kT/b2N) 
ie (20) 
> an aN 
we find 
r=F, eF, (21) 


Because of (18), B > 0. The absolute value of B may, however, 
be much smaller than k7'/b?N . In that case we will find much smaller 
values for B than those found for the exponent in equation (16). The 
time for complete orientation may be as high as Many minutes or 
even hours. 

Thus with plausible values for the physical constants involved, 
we see that the computed time of formation of the mitotic spindle 
and of its disappearance is of the correct order of magnitude. 

If upon division of the chromosomes the substance S is pro- 
duced, which reduces the orienting intermolecular forces, all the fibers 
will contract and the daughter chromosomes will be pulled apart. A 
number of interesting possibilities may be considered, depending on 
whether the f-fibers are first “loosened” from the f,- matrix and con- 
tract alone, or whether the f,-fibers also begin to contract. Consider- 
ing a chromosome, as a very rough approximation as a sphere of ra- 
dius a, — 10+ cm, remembering that it is pulled by an f-fiber which 
has only N/2 links, and considering that the other end of the fiber is 
rigidly attached to the end of the spindle, we find, instead of equation 
(15) 
mila 1. (22) 
b2N 


ar, 
(AON Mab aCe) a 


With the same values of the constants as before this gives us again 
roughly one second for the movement of the daughter chromosomes 
from equator to the pole. This is too short a time. It will, however, 
be considerably lengthened if we consider that the S-substance is not 
likely to abolish the orienting forces at once, but rather weaken them 
gradually. 

One of the quantitative consequences of this theory, which can 
be checked by proper experiments, is the variation of the distance of 
the daughter chromosome from the pole, as given by equation (22). 
Data for that are, however, not available now. 
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The non-transitive character of the peck right relation gives rise 
to different “social structures” in an aggregate of individuals. A meth- 
od is developed for computing the respective probabilities of occurrence 
for each type of structure in small aggregates on the basis of random 
determination of peck right between each pair of individuals. 


Whether sociology can ever become a full-fledged “science” (a 
description of a class of events predictable on the basis of deductions 
from a constant rationale) depends on whether the terms which so- 
ciologists employ to describe events can be analyzed into quantifiable 
observables. 

As long ago as the XVII century attempts were made to con- 
struct a “social physics” (Murchison, 1935). These attempts, inspired, 
no doubt, by the magnificent successes of applying the quantitative 
methods to “natural” (that is, simpler) phenomena, failed, because, 
as Murchison so pithily put it, the “social physicists’”’ made the “mis- 
take of trying to do what physics has done, instead of trying to do as 
physics has done.” The major task of the modern mathematical so- 
ciologist (or the mathematical biclogist, for that matter) is to rectify 
that mistake, that is, to find a way of applying the methods of physi- 
cal science to social phenomena without trying to force them into the 
known frameworks of physics by naive analogies. 

If one asks what physics has done, one might answer that it has 
discovered the inverse square laws of gravitation and electrostatics, 
the principles of conservation of mass and energy and their generali- 
zation, invented descriptions of subatomic events in terms of wave 
and quantum mechanics, etc. But the answer to the question, how 
physics has done, is altogether different. Physics has begun by ob- 
servation and by more or less quantitative descriptions of the results 
of observation. Then it has sought to establish systems of covariables, 
in other words, to reduce the description of events to iteso; then so 
form. Finally, it established a rationale, that is, a mathematical for- 
mulation of “laws” from which all past and many future observations 
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could be deduced. To do as physics has done means to pursue this 
program. As Murchison says, 


The first major task is to identify simple quanta of social phe- 
nomena that by definition are describable in terms of time and space 
categories. After the identification comes the second major task of de- 
veloping a method of measuring these quanta in terms of such cate- 
gories. The third major task involves the application of the technique 
of covariables. The fourth major task has to do with the formulation 
of law and theory. 


Recognizing that the social behavior of non-humans is almost 
certainly simpler than that of humans, Murchison proceeded to iden- 
tify “quanta of social phenomena” in the Gallus domesticus (the 
common chicken), which, as he points out, is a laboratory animal 
with a long and distinguished history. The “quanta” identified by 
Murchison were the so-called Social Reflex No. 1, observed in the 
movement toward another member of the same species, and Social 
Reflex No. 2, fighting a member of the same species. He then estab- 
lished methods for measuring these quanta (Reflex No. 1 by the ve- 
locity of approach and Reflex No. 2 by willingness to fight and the 
success in battle). As covariables Murchison considered mass, mo- 
mentum per second and kinetic energy per second in Reflex No. 1, 
and their relation to success in battle. Another covariable considered 
was “social rank,” that is, the peck order position of the bird in its 
own flock. 


N. Collias (1942) points out that among hens “social rank’ is 
established by the result of the first encounter. He then proceeds to 
analyze the correlations between success in combat with a stranger 
and other measurable characteristics, such as size of comb, degree 
of moulting, mass, and social rank in the home flock. Thus the major 
tasks of identifying, measuring, and correlating “social quanta” are 
under way in the study of events, which may well be considered sim- 
ple social phenomena. In these studies, the unit is taken to be the 
individual, and the observables are his individual characteristics and 
his relations to other individuals within and outside his society 
(flock). 

We wish to attack this problem from a somewhat different point 
of view and by a different method. Our unit will be the society (a 
small flock of birds) and the observables will be the elements of struc- 
ture of that society. Our method will be the method of theoretical 
science: “mental experimentation,” applied to the search for a ra- 
tionale. As has been repeatedly pointed out by N. Rashevsky (1947) 
and others, one need not insist too rigidly on “getting sufficient facts” 
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before attempting to formulate theories, if the “theoretical” nature 
of theories is constantly borne in mind. In fact “mental experimen- 
tation” can often point the way to a fruitful experimental program. 

The search for a rationale is in a way analogous to an experi- 
meéntal program, except that the experiments are performed with pa- 
per and pencil instead of with subjects, materials, and apparatus. 
Even rationales which give “negative results,” ie., fail to account 
for the observed relations, have a certain degree of scientific value. 
Often by observing the “degree of unworkability” of a hypothesis, 
indications are obtained on how to proceed in a search of a better 
rationale. 

In the work of Murchison, Collias, and others, where correlation 
is sought between observables, it is desirable to reduce unknown and 
variable factors to a minimum, since these tend to reduce the sig- 
nificance of correlations. But there is another approach, the prob- 
abilistic, where if the variables are unknown enough, numerous 
enough, a certain randomness may be assumed, and the problem may 
be formulated in terms of random events. We shall adopt this ap- 
‘proach. Since the principles developed here are applicable to any 
analogous situation, we shall speak of societies of “individuals” in- 
stead of chickens. 


General Considerations. 


Peck order or “peck right” is a binary asymmetric relation be- 
tween each pair in a finite aggregate of individuals (a society). The 
relation will be designated by > or by <, so that A > B is read “A 
pecks B” and, of course, A < B is read “A is pecked by B.” The 
relation is asymmetric—if A > B holds, then B > A does not hold, at 
least not simultaneously. Two facts about this relation (as observed 
in bird and other societies) are of interest to the mathematical biolo- 


gist. 

1. The relation is not necessarily transitive—(A > B) - (B > C) 
does not alwaysimply A >C. 

2. The set of relations in a given society may change with time, 
tending to transform the structure of the society into a simple chain 
(Allee, 1931; Murchison, (1935) thus, 


AS BS Ge =e". >) Ze (1) 


The first observation points to indeterminate factors in the dy- 
namics of such a society, while the second points to determinate fac- 
tors. It is, of course, possible to postulate a deterministic dynamics 
and still get intransitive relations. If peck right is completely deter- 
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mined by relations existing between certain inherent or acquired 
characteristics of the individuals, one may formulate the determina- 
tion of peck right thus 


A> R= (AGB). i ie 
BoA BEA) EAs 


where f is a certain function, whose arguments are the respective 
“characteristics” of the individuals involved in the peck right rela- 
tion, and arrows indicate implication. 

Now if f were a function of a single variable, then the met 


cation 


(2) 


A(X) > CYS GY) ee ey rer) (3) 


would have to hold, and peck right would have to be always transitive. 
If, however, f is a function of more than one variable, then it is easy 
enough to construct functions f , such that 


f(A, B) Sf (BGAYE Wey (Bs, Che (CFB ys 
f(C,A) >f(A,C). 


One need only to mark the points (A, B), (B, A), (B, C), 
(C, B), (A, C), and (C, A) ina plane, construct perpendiculars at 
these six points of lengths f(x, y) such that the inequalities (4) hold, 
and pass an arbitrary surface through the six points f(x , y). If 
the equation of that surface is given by z = f(a, y), and if peck 
right between X and Y is determined by whether f(X, Y) = f(Y,X), 
one can see at once that for the individuals A, B, and C the peck 
right will not be transitive. 

One could, therefore, begin the logical development of the peck 
right problem by examining some such function f(x, y) which would 
allow both transitive and. intransitive relations, depending on the 
values of its arguments. The structure of a society will then be de- 
termined by the form of f and by the distribution of the “character- 
istics” among its members. Any change in the structure will have to 
be derived either by taking the characteristics to be functions of time 
or by assuming changes in the form of f with time. 

This approach does not seem attractive, since the number of 
functions one may choose for f is prodigious, and one we no a priori 
reason to prefer one to another. 

The probabilistic hypothesis, on the other hand, supposes that 
the structure of a society is determined to a certain extent by the 
outcomes of:chance events. This means, in the last analysis, that we 


(4) 
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are taking into account our ignorance of the determining factors in 
the dynamics of the society, but we are formulating this ignorance 
mathematically as is done in any probabilistic approach to a problem. 
It is also possible to combine deterministic with probabilistic factors, 
where, for example, one assigns to events probabilities which are 
functions of some known quantities. 

There is another reason besides mathematical convenience for 
introducing chance events into the dynamics of animal sociology. If, 
as is indicated by the observation 2 above, the later stages of a society 
are more “organized” than the early stages, one could postulate the 
initial working of chance and an inherent bias in the situation which 
causes the structure to “gravitate” to a certain form. Thus, if one 
spills peas on an uneven surface, the initial distribution of the peas 
will be nearly random, but eventually they will assume positions along 
the lines of minimum potential energy. 

A natural choice of a chance event is the result of an encounter 
between two individuals. We shall assume that in any such encoun- 
ter, one will enjoy “victory” and the other suffer “defeat.” It will 
be further assumed that the probabilities of the outcomes of such 
encounters can be computed on the basis of the knowledge of the his- 
tory of the individuals and their inherent characteristics. The ‘‘struc- 
ture” of a society, that is, the distribution of social rank will be sup- 
posed to be determined by the outcomes of such encounters. 

One can make several different hypotheses about factors influenc- 
ing the outcomes and about the factors resulting from them. The 
following are examples of such hypotheses, not necessarily consistent 
with each other. 


1. The result of the first encounter between any two individ- 
uals is equiprobable (probability for victory = 1/2 for each). 

2. The probabilities of the results of the first encounter are 
functions of the respective characteristics (independent of time) of 
the individuals concerned. 

3. The probability of the result of each encounter depends on 
the result of the immediately preceding encounter between the same 
individuals. 

4. The result of each encounter depends on the total history of 
encounters experienced by the individuals concerned. 

5. Relative social rank between two individuals is completely 
determined by the result of their first encounter. “hs 

6. The result of each encounter determines the relative social 
rank between two individuals, but other encounters may ensue in 
which the relative ‘social rank may be reversed. The probability of 
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the reversal depends on the difference in the social rank between the 
individuals. 

8. The probability of an encounter taking place at all is a func- 
tion of the relative social] rank. 

Ete. 


One could begin with any self-consistent combination of these 
hypotheses and derive its implications. These implications would be 
conclusions about the probability distributions of various “types” of 
societies, and, if the frequency of encounters were also known or as- 
sumed, one could (in principle) derive expressions for the changes 
of these probability distributions with time. One should thus obtain 
a “wave” of probability distribution. The abscissa of this wave would 
be a particular type of structure of a society (to be defined below) ; 
the ordinate, the probability frequency of this structure; and the 
whole wave would be moving into the third (time) dimension, chang- 
ing its shape as it moved, so that after a long time it would exhibit 
a sharp crest at the type of structure toward which the society gravi- 
tates. 

Thus a surface would be generated by this probability wave which 
could be considered as a complete representation of the “statistical 
dynamics” of the society. 

It remains to define the “type of structure” or simply the “struc- 
ture” of a society. Evidently a definition is desired which would en- 
able us to recognize structure by observation. In a society of N indi- 
viduals, each will have N — 1 peck right relations of which r will be 
dominant (the right to peck) and N — 1 — r submissive (the neces- 
sity to be pecked). There will then be a distribution of numbers 
(7,, %2,+*++ fy) among the N individuals such that 


N 

Dri = I1N(N—1). (5) 
This set of numbers (7,, 72, +--+: ry) together with all its permutations 
will be defined as a structure of the society. A structure can also be 
represented by a diagram where each individual is placed on a level 
determined by the number 7; of dominant peck right relations he en- 
joys. Thus a society of four individuals can have exactly four struc- 
tures, diagramatically represented in Figure 1. Arrows indicate peck 
right; the numbers indicate the dominant peck right relations of each 
individual. The problem of determining the dynamics of such a so- 
ciety, given a set of assumptions about the occurrence and results of 
encounters, reduces to the calculation of the probability of occurrence 
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FIGURE 1 


of each structure at each time in the history of the society. In this 
paper the following simple assumptions will be made. 


1. The results of the first encounter between any two individ- 
uals are equiprobable. 

2. The result of the first encounter determines permanently the 
peck right between the individuals concerned (i.e., the winner will 
always peck the loser). 


Thus the shape of the probability distribution wave is perma- 
nently established (independent of time) as soon as all the 4N(N — 1) 
encounters have taken place. An alternative interpretation is that 
we are limiting our observations to the period when the structure 
depends on the results of first encounters only. 


Three Individuals. 


Since only one type of structure is possible for a society of two 
individuals, namely, (1, 0), we shall begin with the case of three 
individuals. 

There are 2? = 8 possible sets of outcomes of the three encoun- 
ters. Each outcome leads to one of two possible structures, (2, 1, 0) 
and (1, 1,1). The respective probabilities of these structures are 
3/4 and 1/4, since six of the eight outcomes map on the first (sim- 
ple chain), and two on the second (simple cycle). This solves the 
problem of the three individuals. 


Some Aspects of the N-Individual Problem. 


One would like to generalize the solution to N individuals. How- 
ever, difficulties are encountered even under the extremely simple 


assumptions made above. 
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The probability of an event is defined as the ratio of the number 
of ways the event can happen as a result of certain occurrences to the 
total number of possible occurrences, assumed equally likely. Now in a 
society of N individuals there will eventually take place 4N(N — 1) 
first encounters. Since there are two possible results to each encoun- 
ter, the number of all possible occurrences (where an occurrence is 
the complete set of results) will be 2?%*. To compute the prob- 
ability of a certain structure, one must calculate how many of those 
“sets of results” will “map” on a particular distribution of peck 
rights (71, 72, °°: fy) or on any permutation thereof. This would 
be a simple matter if there were a one-to-one correspondence between 
each set of results and each ordered set (11, T2, -**+ ?y). Then the 
probability of the structure (7, , 72 ,--:- 7y) would be simply the num- 
ber of distinguishable permutations of the set divided by 2?%?-», 
That the mapping is not in general one-to-one is shown by the fol- 
lowing counterexample. 

Each set of results can be represented by a skew-symmetric ma- 
trix, for example, in a society of four individuals by 


7 ig aa 3 Cent 37) 
Aj} 0 -1 -1 41 
B | +1 0 +1 —-1 
c|/+1 —-1 0 —-1 (6) 
D|/—1 +1 #41 0 


where the article a +1 in the X-row, Y-column indicates X > Y and 
a—l, X < Y. Note that the matrix (6) maps on the structure 
(1,2,1,2). But so does the matrix 


sAle Reset Cala) 

Pte a ee 

BA gels Ae hee 

Peg Raat te oer RET Pe a (7) 
D.|| ct lanersal Agee a LOF 


Without solving the problem of determining the probability dis- 
tribution of structures for all N, it is nevertheless possible to re- 
word it in such a way as to state the results in terms of properties of 
such skew-symmetric matrices. In this form the problem may, per- 
haps, prove susceptible to attack by methods derived from the theory 
of these matrices. For this reason, we introduce an alternative defini- 
tion of structure. Let the structure of a society be defined by a set 
of numbers 


(11' 5 Te! +++ Ty’), | (8) 
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each of which represents the excess of the number of dominant rela- 
tions over the number of submissive relations associated with an in- 
dividual. Then evidently 


N 
= r; =0. (9) 

Theorem I. In a society of N individuals the structures defined by 
(ri, Yo, +--+ ry) And (1ry', Ye’, ---- ry’) are in one-to-one correspond- 
ence. 

Proof. Let a given structure be given by (1,72, «:* ry) in the 
sense of the first definition and by (71', 72’, --:- 7y’) in the sense of 
of the second. Then 

rr, —?T,— (N —1—7;) —27,—N +1. (10) 
Also 
r,—1/2(r; + N—1). (11) 


The right-hand side of equation (11) always represents an in- 
teger since from (10) we see that 7,’ is even or odd depending on 
whether N is odd or even. Thus the 7,’ are uniquely determined by 
the 7; and vice versa. This proves the theorem. 

Consider now the class ™M of all skew-symmetric matrices 
(a:;) of order N, all of whose non-diagonal elements are units, that 
-is, d;;—= 1 or —1; 4;, = — 4;;. 

Theorem 2. Let (ai;) be a matrix of the class M, and let 


j=1 j=1 


N N N 
S.=( D A+ Ds 423, o Baw] 


be the ordered set of its row sums. Let p, be the number of distin- 

guishable permutations of the set S,, and m, the number of matrices 

in M giving rise to the set of row sums S,. Then the probability of 
N 


the structure (r1', Yo’, ---- tx’) , where ri’ = > ai; ts given by 
4 j=1 


P(r , 12! 5+7°* Tx) = 29 Pym « (12) 


Proof. Each matrix (a;;) represents an equiprobable set of en- 
counter outcomes. Furthermore, the row sums of the matrix (aij) 
represent the excess of victories over defeats of each individual and 
therefore, according to our assumptions, the excess of his dominant 
over his submissive relations. Therefore m; represents the number 
of sets of outcomes mapping on a particular ordered set S,. But any 
permutation of S;,, S;‘ can be interpreted as simply a re-labeling 
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of the individuals. Hence the number of matrices mapping on 
S, (G= 1, ---- py) is likewise m,. Moreover sets of outcomes map- 
ping on distinct sets S are certainly distinct. The total number of 
matrices in M for a given N is 2!%", The number of distinct ma- 
trices mapping on the wnordered sel S;, is pxm, . This proves the theo- 
rem. 

Because of the one-to-one correspondence between the structures 


(71, T,-°-* Ty) and (71, %2', -+++ fy’) established in Theorem 1, we 
have the 
Corollary. The probability of a structure (11, Y2, +--+: Yn) ts also 


Q3N(N-1) Diy, « 

A general method of finding m, for any S;, and any N is at this 
time unknown to the author. However, for N not too large, say 
N & 6, the probabilities of all structures can be easily computed by 
a systematized procedure, which will now be described. 

First let us consider a procedure for writing down all possible 
structures (71,72 ,°*:- 7'y) of a society of N individuals. The simplest 
such structure is a simple chain (complete hierarchy). It is repre- 
sented by the set of numbers 


(NON 2 she ote care 0). (13) 
This structure is a mapping of the matrix 
OT eee eee 
pee Ue es 6 ae 
rt) + 
(14) 
—— a eee eC) 


where all the pluses are above the diagonal (the units have been re- . 
placed by their signs for simplicity). All other matrices of class M 
are obtained by changing a set of pluses in matrix (14) to minuses 
and the symmetrically situated minuses to pluses, i.e., by a set of in- 
terchanges of pluses and minuses. But not all such fierciioes will 
give rise to structures different from (13), because some of them may 
map on permutations of (13), which, by definition, represent the same 
structure. Such permutations can be avoided if we consider uy, such 


interchanges which preserve the order of magnitudes in the > j=", 


n=1 


that is, their non-increasing order. 
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Such interchanges can be represented as follows. Consider N 
boxes containing respectively N — 1 ,N —2,---- 0 balls. We are in- 
terested in the operations which throw the set N — 1 ,N —2,-:---0 
into another set (71, 7, ---- 7vy), where Yr, 2 r; fori < 7, and such 


N 
that S| r; = 4N(N — 1). Such operations can be accomplished by 


421 
transferring the balls in several ways, but if the balls are not dis- 
tinguishable, for each transformation 


Ni Nat 0)? (715 Te, 2x) (15) 


there exists a unique, most direct way of transferring the balls, where 
each box either gives up balls or receives them or neither, but not 
both. For instance the transformation (3,2,1,0) > (2,2,1, 1) 
can be accomplished by transferring a ball from the first box to the 
second, from the second to the third, and from the third to the fourth. 
But there is a unique direct transfer with the same result, namely 
where a ball is taken from the first box and put into the fourth. It 
is also evident that all such direct order-preserving transfers are 
from left to right. But each transfer represents an interchange of 
signs in matrix (14) and also avoids the permutations of (7;). Hence 
we have 

Theorem 3. The number of possible structures of a society of N 
individuals equals the number of ways direct transfers of balls can 
be made preserving non-increasing order in a set of N boxes contain- 
ing N — j balls respectively (j =1.----N). 

By systematizing these transfers in some convenient manner, we 
can write down all possible structures of societies of N individuals, 
if N is not too large. Thus a society of six individuals gives rise to 20 
possible structures. Beginning with the complete hierarchy (5,4, 3, 
2,1, 0) the other structures can be written down by imagining all 
transfers of one ball, then all transfers of two balls, etc., as shown 


in Table 1. 


TABLE 1 
S,: (5, 4, 8,4, 1,0) S,.2°(5, 4, 891, 1, 1) 
So: (4, 4, 4, 2, 1, 0) Sic: (4, 3, 8, 8, 2, 0) 
Si (4,4,8, 5/110) S eii(4; Sy onesds 1) 
S,: (4, 4, 3) 2, 2, 0) S(4,5)952,2) 1) 
Se (Aj4,8y2, 1,1) Sus (4d 2 251) 
S,: (5,8, 8, 3, 1, 0) S.g2 (5; 8, 2 2, 2, 1) 
Se (Bro, 8) 2,260) Siz: (3, 8, 8, 8, 2,1) 
S,: (5, 8, 8, 2, 1,1) 5,308, 8, 8,2, 2, 2) 
S(64;2,.2; 4,1) Syoi (5, 2, 2,2, 2, 2) 
Bild, 452,21, 1) Se, xl 4p Sy thoy 22) 
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It remains to indicate a method for computing the probabilities of 
structures. We have seen that these probabilities are determined by 


the number of matrices of class M which map on the set (71, 72 ,°°:: Tw) 
or any of its permutations. Really it is necessary only to determine 
the number of matrices mapping on the ordered set (71, 72, -*:- Tw), 


as indicated in the proof of Theorem 2. The total number of matrices 
mapping on a structure will then be the number mapping on an or- 
dered set (71, 12, ++: Py) times the number of distinguishable per- 
mutations of this set. 

The society of three individuals has already been described. For 
the case of four individuals, we have 2° = 64 possible matrices. The 
ordered set (3, 2,1, 0) can be obtained only from the matrix 


Oe EER ea 
— 0 


| o+ 


+ 
ee 
—- = 0 


This set has 24 distinguishable permutations. Hence the probability 
of the complete hierarchy in a society of four individuals is 


P(3,2,1,0) =24/64=38/8. 


The ordered set (2, 2, 2, 0) can be obtained in two different 
ways, 


(Pate eee ee 0}, ticanieny shit bck 
cmasicige Ocontemet eines sink 20) tote oiheete 
eee car. wee eee —- + 0 + 


Its distinguishable permutations are four. Therefore 
P(2,2,2,0)™=2X4/64=1/8. 
The set (2, 2, 1, 1) can be obtained in four different ways 


0 = te _— 1) + at _ 
ee A AONE Ra —W0.a he oe 
re ae 0 Ae i+ —_ 0 = 
pw ee tes SO Se ee 
0 = -- cus 0 — + ae 
ats 0 aie — + 0 a. ati 
Berto aX 0 se — + 0 = 
2.2.58 Aer ee oe ee 


ASATOL BAPSFORT 1% 
P(,2,1,1) —4*6/aA-=3/s. 
Vimally (31,1, 1) (an be difzined in two wars 


——s SS ~ Vo = 


a oe a a 
ee ee ts 
i — % * | —— + of] =a 
_ + s 6 | he — - 0 | 


PG ,.1,1,1)=—2* 4/A-—1/8- 


/ This completes the case of four individeels. The probchilitics of 
the structures occurring in socicties of five individuals con be cosy 
Comspuied by this method. The results are given in Table 2 


TAZLZ 2 

® Meret Wiser Noeezz oF a 
Srzeacrizz os Orsezza Ger PremerreTios Pepezemzrr 
i 12 22x 

2 v4 227 @ 

4 BB 25> 13 

z B 2* ~ 23 

2 BB 22 > #B 

Z B 28 523 

A 2 2-2 x & 

1“ wy 2 > Zs 

a a ee ae 
The foregoing results muplicitiy predict the Gtrowtis of Sree 
wes in 2 lerge member of socicties Af 2 given wumber of individuals, 

the hypotheses used im their derivations are vad On the 

her hand there are strong indécctioms thet in any actmal society Tike 


— imam : 
LL The hove ressits can be used 28 2 theorctacal Seppe ome 

‘ oe ere: noes me 7 
; ese ee 

= a? 6G : _— 52) 5 ae : : 
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mined by the result of the first encounter, the probability distribu- 
tion will depart from that predicted here so as to weigh more heavily 
those structures which are closer to complete hierarchy. For exam- 
ple on the basis of random results of encounters, the probabiliy of a 
complete hierarchy in a society of three individuals is 3/4. With the 
introduction of some bias, we would expect it to exceed that value. 

3. Whereas the statistical work of correlating physical char- 
acteristics and social rank with success in combat is most significant 
when the opponents are not evenly matched, the results of the pres- 
ent theory can be best realized when they are matched as evenly as 
possible. Thus experimental work can be extended to another range. 


This work was aided by a grant from the Dr. Wallace C. and 
Clara A. Abbott Memorial Fund of The University of Chicago. 
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CONTRIBUTION TO THE LINKAGE THEORY OF 
AUTOPOLYPLOIDS: II 


HILDA GEIRINGER 
WHEATON COLLEGE 


In the first part of this paper, hereinafter referred to as I, a general 
segregation distribution was introduced for autopolyploids (2s-ploids) 
with m loci, and r alleles. Random mating, chromosome segregation, 
distinct generations, and equal segregation distributions for males and 
females were assumed. In this second part, using this segregation dis- 
tribution as a basis, a recurrence formula is established which enables 
us to compute the distribution of gametes for any generation, if this 
distribution is known for an initial generation. This initial distribu- 
tion of gametes is derived from the initial distribution of genotypes. The 
limit behavior of these distributions is completely described in a gen- 
eral limit theorem which contains as particular cases the limit theorems 
for diploids with m loci, and for 2s-ploids with one locus (Geiringer, 
1944, 1945, 1948). 


1. The General Recurrence Formula. 


The general recurrence laws for polyploids with m loci will now 
be derived. We write, as before, for the distribution of gametes and 
for the distribution of zygotes in the mth generation p™ = p, 
w™ =w, and p') = yp’, w™) = w’'. We consider a 2s-ploid with 
m loci and 7 alleles. We denote the ith allele corresponding to the kth 
locus by ay, (i =1,---- m, k =1,---- 7) and the distribution of 
gametes by 


p (a,,* a elere C774" pote Qina?™ nes Omr?™*) : 
Zig tee + 2ir=8 
= ik rp OQsiG Nd 
If r > s some of the 2; will necessarily be zero. There are 
s! s! 


Bal<nseae! Zt = sekmes 


such gametic probabilities for a set of zx and we want to express 
S = Dl (dart oe Oyy2h 88+ Omi? +++ Une") in terms of the w’ as was 


done in the case s=m=2 (section 5 of I). 
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We shall need the marginal distributions of the p-distribution. 
The definition is analogous to that in the particular case considered in 
section 5 of I. A general example of marginal summation will suffice. 
Consider m=3 loci, r—=4 alleles and s=6. Denote, in this example, the 
alleles by a1, °+:- , 4, 0:,-°-*- , Di, C1 » °°» Ca. The p-distribution 
is a function of eighteen variables each of which takes on four dis- 
crete values. From this distribution we may derive marginal dis- 
tributions of order 17, 16, ---- , 1. There are three of order 17, 
six of order 16,---- , three of order 1. As an example consider the 
marginal distribution of order 11: p(a@,a.b,b;?¢,¢,¢,2). We will use 
summation letters: a taking on the values .a,,---- , a, 6 taking on 
the values b,, ----, b,, and y taking on the values ¢,,---- , ¢,. The 
marginal distribution in question is then defined as the sum of 4" 
terms: 

ey SY DPI P (Q12020, 0203003781826, CoC.? yi 72) 
a1 G2 ay Bi Be Yr Ye 
: "=p (aad, bs?C1C2€,?) . 


It is obvious that in the same way any marginal distribution for any 
m,7, 8 may be defined. 

We return to the expression S + p’ (@:*" , --++ , Gmr*™") considered 
in the first paragraph of this section, which we want to express in 
terms of products of the values of w’ and of thes.d. According to the 
meaning of the three distributions involved, considering our hypoth- 
eses of random mating and of chromosome segregation, each w’ will 
be multiplied by a sum of certain J-values; or, the other way around, 
each A-probability will be multiplied by a sum of certain w’-values. 


We want to find the coefficient of A(e,, ---- , em) in this expansion 
(see section 4 of I for the definition of A(e,,----, €m)) ie., we want to 
collect all those w’-values which multipy a fixed A(e, , ---- , em) SO as 


to have a formula 
Sp’ (4,7 awe Gre ponte Qmi?™ ones Omr?™") 


=O (epee tee lonne be (eS 

Here, according to the meaning of the s.d., a A(e, yrs €m) Is mul- 
tiplied by a sum of w’-values each of which contains before the semi- 
colon «; “correct values” corresponding to the ith locus and after the 
semicolon (s—e;) “correct values” corresponding to this locus. By ¢&; 
“correct values” we mean ¢; of these values which stand with respect 
to the 7th locus on the left-hand side of the preceding equation which 
means ¢; values of the set a;,°" ---- ai,”* (( =1,---- m). Hence such 
a w’ will be of the form 
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Uw (a, 44 Q4,7" He Qmnx™™ ny? ™? ‘ 
444 M7 9 Ayn! 2+ Qn! ™* ): 
where 
Lig e+ Lip =8;, Yir'F +9 Yip = 8 — 853 
Lik T Yin = Zix » Zig + +++ Zi, = 8, (2) 
Gis; em), (k= 1,-++-7). 
In such a w’ there are still (s — e,) + (8 — &) +--+» + (8s — em) 
“free” or “empty” places before, and e¢, + & + - wath Em “free” or 
“empty” places after, the semicolon. Consider for the moment, to fix 
ideas, a fixed choice for the values on those LS nes ee By Jose sees 


[(s — €m) + en] = ms empty places. Let, e.g., the first s of them all 
be equal to a,,, the second group all be equal to a», , etc., and let us 
count the number of w’ values which multiply 4(e, , ---- em) for this 
particular choice and for a particular set of xi, yix. The e; “correct” 
values before the semicolon can be arranged in( . 


a 


)ways on their s 


a ways 
See y v 


on the s places available for them. Moreover, these «; values can 


places and the (s — e;) values after the semicolon in ( 


e;! ee. 
still be arranged in ———_——— ways within the chosen «; places and, 
Vial --++ Li! 
se (s— 2) 1 geet, , 
similarly, there are SSE rot possibilities for the (s — «;) “cor- 
ins 7o** Yir ! 


rect” values after the semicolon. Hence the coefficient of a w’ is 


5 Ree s 2 e;! (s— &) ! 
i ! a rise ‘i 
i &j Lizsrres irs Yire ees Yir- 


N ex, according to equation (4) of I replace each w’ by the cor- 
responding product of two p-values and divide both sides of (1) by 
“TT 1 


Cot) IT eer The coefficient of 4(€,, --:* €m) in the expan- 
sion ge - (117% ++ ++ Gn,*"") is then a sum of terms of the form 
1 ee ( 8 )( s |) eter tl | 
(s!)™ i &i bp] ty) - Liz! Yal-<s<Yir! 4 6) 
D (yy +++ Oye + 5 Osyep Bat 9 Gig g OF == 5 «) 
BO (POLO | aa Dyn l™1 228 * Amr! ores). 


So far we have assumed one fixed choice of values on the (m:s) 
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“empty” places. Now we abandon this assumption. In the first 
p-probability in (3) there are (s—e,) empty places after a,"",----, 
and (s—e,) empty places after @m,?"*. Similarly, there are «, 
(i =1,---- m) empty places in the second p-probability. On the 
(s—e;) free places in the first p-function we have to write the r values 
Qix, ««-+ Gi, in all possible 7*-* combinations, and, similarly, for the 
second p-function. Then all these products of p-values are to be 
added. This amounts to an [(s — e:) + e:] +--+: + [(s — em) + Em] 
= ms-tuple sum of such products. But this sum of products obvi- 
ously resolves into the product of two sums, one an [(s — ¢,) + -:::- 
+ (gs — e,) ]-tuple sum, the other an [e, + ----+ €,]-tuple sum, these 
sums being nothing else than the marginal distributions 


DAA eae s Bir os ones @m1 ..0e z ) 
D(A. Qir*s Ama ™ Qmr ™" 


and 
/4) (4,4 euae Qy pir POSS Qin t™ eeee Qnr’™) 


of order [e, + ---- + em].and [(s — e,) + ---- + (s — em)] respec- 
tively. Thus, for a fixed set of xix, yix the coefficient of A(e,, --+- Em) 
becomes 


n ( =i 4 ej! ; (3 2;)1 Bighees: Sir! 
i gj Gj / Hiylsss+Liel Yl Yir! s! 


“Pp (4,7 ses Amr ine) 24d) (43 *9* Onr ine) ’ 


(4) 


with (2) holding for the Vik * Yik » Zik» 
This last expression can be simplified. Observing that yi, = zix 
— “i, we have 


s! 6;! (S—6;) 1) 2p lee! 
€;!(s—e;) ! Pr lene tee (tet otey ate s! 

ae ae Zirt 

 fal@n—@n)! al (rte) 


Hence, introducing as before (section 4 of I) 


Tae =| :) he ( o) Alena em), 


pa a! 
 Bla— py 


we arrive at the very simple and general result 


and using as always the notation ( : ) (0OSP=a); 
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2... z 2 2 
(diz OO 8+ Chany eee i= Lany =) 


== » Soi a ley: -- €;,) 


(5) 
BeBe Bee (7 4 fir ae emi BAG &mr 
T34 yy Tmt tmr was v2 f Uma i 
0) (a437 cee Qmy™™") -p (G0 oidass Qing? ™ 7 r) : 
where 
fuer **-2i,— 8, Hohn, SOOM yea ae 
(6) 


OS ti ae ((=1,----m), (kK=1,----7). 


This is the general recurrence formula for the distributions of gam- 
etes in two successive generations for a 2s-ploid with m loci and r 
alleles under random mating and chromosome segregation. It con- 
tains as particular cases the respective recurrence formulas of Geirin- 
ger (1944) and (1948a). 

We finish this section by considering a few particular cases. 
First, let us consider a homozygous organism. We then put 


Qyy 08 = Oy = Ag ye ny 8 Amr =An, 


and obtain the simple result 


0....¢ s 0 


acon s 
p'(A,o+-+- An’) = = oe SS Ue, ) -( ).- 


Em 


= p(A,* cae ee Ane”) ae (A,3-*: Ca Am ™).s 


The interpretation of (7) is obvious: [(e,-+-- én) - ( . ) tees Es) is 
the probability that «; maternal and (s — «;) paternal genes will be 
transmitted with respect to the ith locus (4 = 1, ---- m), while 
(p(A,® «++» An") is the probability of a maternal gamete with «; Ai- 
values, and p(A,*-® ---- A,,-&™) is the probability of a paternal one 
with (s — e;) A;-values, in the nth generation. Under random mat- 
ing, the sum over all ¢,, ---- &m of the products of these three prob- 
abilities is the probability of a gamete of type (A,' ---: Am’) in the 
(n+1)st generation. It is not difficult to give an analogous inter- 
pretation for the general result contained in (6) (see Geiringer, 
1948a, page 259). The biological situation considered in the present 
paper (general m,s, and 7) is in a certain respect the most general 
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one in which the recurrence formulas have the clear and simple struc- 
ture exhibited in (6) or (7). If chromatid segregation is studied we 
find recurrence laws of an essentially different, much more compli- 
cated type, as will be shown by the author in a forthcoming paper. 

Next consider the important case of just two alleles, r=2 at 
each locus, A; and a; (i=1,-:-- m). Formula (6) then yields the 
result 


p’ (A,710,,9-4 siaels Am?" Qm>*™ ) 


= FS Ue: +: em) 


~ &1 SS L4 Tue em S—<a 3 
pe ON ae alee (8) 
-p(A,™a,2™ sees Aye Ame ™*™ ) 


-? (A, 21g 2 Fe ers eons Ay?" Dy S2m-Em +m) ‘ 


In (8) the summation with respect to the x’s is reduced to an m-tuple 
one rather than an 2m-tuple one, as would follow from (6). 

Finally, and this will give a very useful result, consider the case 
where r 2 s and 


Say tt yy S Boy SH Sop SH SS St Sm = 1. 


Without loss of generality we obtain the formula 


0....8 


0....8 
D' (Arya +++ Ag 2+ ** Amy ++ Ams) =>. Sh E(e,, +--+ &m) 
om 


1 


- [Dp (aQi +++: @ 2886 Amie @ ) (9) 


1&1 MiEm 


*P(Are stig s+ Ome 4a -2e* Ome) 3 oe 


There are in each bracket ( Jon ( “ ) products of which the 
1 m 


first has been written down. Let us however mention that the bracket 


which multiplies I(u, uw, ---- w) contains 4. ( : Me different terms 


only, while that whole bracket is multiplied by “two” [see (13)]. 
Let us show an illustration of (9). Consider m=2 , s = 8, write 

@,, A, ds for the first three and b,, b., b; for the other alleles. We 

find ; 
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D' (4,d203b1b2bs) = [1(00) + 1(33)] p (a,4,5b,b2b,) 


+ [1(10) + 1(28)] - [p (a) p(cuagb,b.b,) tee hoe] 
+ [1(01) + 1(32)] - [p(b,) p(a,a,a5b.b,) +--+ .--] 
+ [2(20) + 1(13)] : [Pp (G12) p (asbibob;) oe Sea re ond | 
Ta (hl) 1(22)].- [Pp (a1b1) p (a.asb2bs) (10) 


+ +--+ (nine terms) ] 
+ [1(02) + 1(31)] - [p(b1b2) p(Qidedsbs) +--+- fee] 
+ [1(30) + 1(03)] p(a,a2a5) p (b1b2bs) 
+ 1(12) + 1(21)] - [p (a...) p (a2a3b3) 

+ +--+ (nine terms) ]. 


Next consider the case m=2 , s=4. We want to write down only the 
term which multiplies 2/(22). It is: 


21 (22) [p (4,a2b,b,) p (a30,b3b,) + p (@,a2b,b3) 
-P(d;0.b2b,) + p(G,a2b,b,) p (a,a,b2b3) (11) 
+ .---- (eighteen terms) ]. 


Since [(¢,, ---- &m) = U(s — &,, ---- $ — &m) the formulas for 
s—=2u as well as for s=2u + 1 take on the following form 
Oreo fe 
p (Ay see Aig tt * Ami tt*? nt) =<? oe l(e, AiO Em) 
' [DP (Ga +++ + Qe +++ Ami *+** Ome, ) (12) 


“DP (Gre 41 ses Aig tees Am,e, + ne es Ging) = v2] 3 


Here >’ indicates that for s=2u the expression in the brackets which 


’ 2u\|" a 
multiplies 2/(u, w, ---- w)contains $ - ( iu )| terms only. As an 
example consider for m= , s=2: 


D' (A,02b,b56,C.) = 21(000) p (a,d2b1b26;C2) 


+ 21(100) [p (41) 7 (@2bib2C,02) + -++-] eves be 
+ 21(110) [9 (@,b1) p (deb26102) eves bees Hee] (13) 
HS a Bt 


+ 21(111) [p (arbic,) p (deb2C2) eves teers Heese]. 


From (9) or (12) we can now find any formula regarding partly 
equal alleles by equating some of the a, to each other, some of the 
a., to each other, etc. If for example in (11), a:=a.—a;—a,—0 the 
coefficient of 1(22) will be: 


12 1(22) - [p(a?b,b2) p (a?b,b,) + p(a*b,bs) p(a*bebs): 
+ p(a*b,b,)p (a?b2bz)] . 
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It is thus easy to write down any particular recurrence formula. By 
means of such a formula the distributions of gametes (and therefore 
because of equation (4) of I the distributions of zygotes) can be 
found for any generation if the initial distribution is known. 

The recurrence formulas are, mathematically, quadratic differ- 
ence equations whose constant coefficients are the values of our s.d. 
The decisive role which the s.d. plays in the organization of our 
biologically meaningful recurrence relations seems to hint that our 
concept of a s.d. is not inadequate biologically. On the other hand, 
once our recurrence formulas have been established, they are com- 
pletely independent of any particular assumptions regarding the para- 
meters in the s.d. Random chromosome segregation or many other 
models of segregation may be used (see also section 4 of I). No re- 
currence formulas have been given by other authors concerned with 
our problem. 


2. Recurrence Formulas for the Marginal Distributions. 


If we wish to integrate the recurrence formulas of the preced- 
ing section or if we want to study the limit behavior of our distribu- 
tions, we need recurrence formulas for the marginal distributions. 
They can be derived from the recurrence relations of section 1 of II 
(see also section 5 of I). Consider first the case m=s=r—=2 where: 


p' (a,0,b,b,) = 21(00) p (a,d2b,b.) + 21(01) 
-[p(b:) Pp (Qyaeb2) + P(b2) DP (a1d2b,)] + 21(10) 
é [p (a1) P (a2b,b2) ay Dp (a2) p (a,b,b2) ] 
+ 21(02)p (d,d2)p(bib2) + 21(11) [p (arbi) p (dab) 
+ p(a,b2) p (a2b,) J, 
and 


p' (a,0.b,?) = 21 (00) p (a,@2b,7) + 41 (01) p(b:)p (@,@2b,) 
+ 21(10) [p(a.) p (aeb,2) + p (a2) p(a1b,") J 
+ 21(02) p (02) (b,?) + 4l(11)p (a,b,)p (a2b,). 
Upon addition we obtain 
p' (a,020;) = ' (a,a,b,b2) sah fs (a,a,b,b,) 
= [21(00) + 21(01) ]p (a,a2b,) 
+ ['21(01) + 21(02) ]p (a,de) p (61) 
+ [21(10) + 21(11)] [p(a,)p(ab,) + Dp (a2) p(a,b,) ]. 


Here the coefficients of p(a,a.b:) or of p(a.a2) p(b,) to the right are, 
of course, the marginal distributions of our s.d. as introduced in sec- 
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tion 4 of I. Using the notation of section 4 of I we have to write 
D' (41d2b1) = 21 (00) p (ay4ab1) + 21 (01) p (ads) p (0) 
saan Oye) (10) [p (@,)p (a2b,) + p (@2)p (a,0,) J, 
with the symmetry properties as in equation (36) of I 


LG) (90) =1M (21), 1 (10) =lo~™ (41), 
§69) @) (01) — [0 G) (20) - 


(14) 


We may derive a recurrence formula for '(a,b,b.) in the same way. 
N ext, by a twofold summation, still for s=m=2 we obtain 


p' (d,d2) = 21 (00) p (a,4) + 20 (10) p(a.)p (ae), (15) 


where 
12) (00) =71(00) + 21(01) + 2(02) ; 
Boe 10 y= 1 (10) = 27011), $12), 
and 


p' (a,b,) = 21 ™ (00) p(a.b,) + 219 (10) p(a1)p(b1), (16) 
with 
14) (00) =1(00) +12(10) +1(01) + 1(11); 
1G) (10) =1(10) + 2(20) + 1(11) + 1(21) ; 


and, of course, in a similar manner we obtain the recurrence rela- 
tions for p’(b,b2), for p’(a,), and for p'(b,). 

It does not seem necessary to show the marginal character of the 
distribution to the left in an explicit way in the notation, since this 
becomes quite clear by the superscripts in the values of the s.d. to 
the right. If, for example, the distribution to the left is p' (a,b1) 
while on the right there is a term containing /) (00), it is clear 
that we deal with an s=m=2-organism. Or, if in a formula where 
p' (a,a,b,b.c,) stands to the left, on the right there is a term contain- 
ing, for example, 1 © (0000) we see at once that the organism 
considered is of type m=4 , s=5 since 2+3=2+3=114—015=5. 

Using the concepts introduced in section 4 of I and the resuits 
of section 1 of II we obtain easily the recurrence formula for any 
marginal distribution. Consider the case just mentioned and write 


now for the sake of brevity 1“ © = 1. We have 
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p' (A,%b,b2¢,) = 21 (0000) p (4102b,05¢;) 
+ 21(1000) - [p (a1) p (a2b102¢,) + DP (a2) p (a10,52¢;) J 
+ 21(0100) - [p (01) p (Gia2b2e1) + p (bz) p (a142b;¢,) ] 
+ 21(0010) p (c1) p (a,02b,b2) (17) 
+ 21(2000) p (a,d2) p(b,b2¢:) + 21(0200) p (bibs) p (Aid2¢1). 


If in (17) we then put a,=a, we find a formula for p’(a,°b.b2¢1). 


A general formula for the marginal distributions, analogous to 
(12), is as follows. Put s—»i=2u;:, or s—i—2ui+1 respectively 
(i=1,---- m). We obtain 


p (G4, Pode Qi,s-v, seiate lone owes Gm,s—v ) 
Ow... pa 0...6pLm 
SO ace Sar es ey.) 
E4 Em 
(18) 
. {p (a, ahead Qi,2, SOOKE ROUT Ome) 
-p (die, cs Do a Os sy aS Om, pate Pay On, sv ) sa gekhed -] - 


1 


The recurrence formulas (12)—(18) are quite general and seem 
to reveal a regularity of the heredity mechanism. They are so simple 
that we may write them down for any m, 7, and s without computa- 
tion. The values of the s.d., or of the marginal distributions of the 
s.d., appear each time as the constant coefficient in these quadratic 
recurrence relations (finite difference equations). The essential con- 
cepts are: (1) thes.d., and (2) its marginal distributions; (3), the 
d.ga., and (4) its marginal distributions. Without the use of these 
concepts the computations become hopelessly involved, as exempli- 
fied in the explicit computation of section 5 of I. Recurrence for- 
mulas for the simplest cases are due to Jennings for s=1, m=r=2, 
and to Haldane for m=1 , s=r=2. The essential restrictions for our 
recurrence laws are: (1) random mating, (2) chromosome segrega- 
tion rather than chromatid segregation, and (3) the s.d. independent 
of the genotype of the organism and of “n.’’ (It may be noted that 
the essential structure of the recurrence laws does not change if dif- 
ferent s.d.’s are assumed for males and for females.) 

It is not difficult to work out an integration of our recurrence 
relations and to arrive at explicit formulas in a few simple cases. By 
that we mean, for example, that p (a,a.b,b.) is expressed in terms 
of p (a,a.b,b.), of the marginal distributions of p™, of the para- 


HILDA GEIRINGER 207 


meters introduced in the definition of the s.d., and of course, of “n.” 
The explicit formulas for m=2, s=1 and s=2 , m=1 are due to Jen- 
nings and Haldane, and for s=1 , m=3 to the author, who worked out 
a general procedure for s=1 , any m, and for m=1, any s, (see Geir- 
inger, 1944, section 8; 1948a, section 6). 


3. The Limit Theorem for m=s—2. 


We begin with the consideration of the case m=s=2 where all 
possibilities can be discussed explicitly. This discussion will be need- 
ed in the treatment of the general case. 


1) With the notations of section 2 of I we first assume 


1,(0) =2(00) + 21(01) + 1(02) <4; (19) 
1,(0) =1(00) + 21(10) + 1(20) <4. (20) 

It follows that 
l.2(00) =1(00) +7(10) +2(01) +7(11) <4. (21) 


In fact, if 21,.(00) = 21(00) + 21(10) + 21(01) + 21(11) =1 would 
hold together with 


21(00) + 42(01) + 41(10) + 41(11) + 21(20) =1, (22) 
it would follow that 1(10) = 1(01) = 1(20) =1(11) = 0 and 1(00) 
= 1/2, and then (19) and (20) could not be true. We formulate this 
intermediate result explicitly: 

7 00) 31/290 —1,6(10) It. 1 (01) =U , then 1,(0) =1.(0) 
= 1/2. Hence if at least one of the two inequalities (19)—(20) holds, 
Cio > 0 . did: igs 

Now consider the recurrence formulas. Writing p” = p, 
p(n+1) gy. CRAP kth abe bm, == OE 

p' (AaBb) = 21(00)p(AaBb) + 21(10) - [p(A)p (aBb) 
+ p(a)p(ABb)] + 21(01) [p(B)p(AaB) + p(b)p(AaB)] 


23) 
+ 21(20)p(Aa)p(Bb) 
+ 21(11) [p(AB) p(ab) + p(Ab)p(aB)], 
and 
p' (AaB) = 21 (00) p (AaB) + 21 (01) p (Aa) p(B) a 


+ 21 (10) [p(A)p(aB) + p(a)p(AB)]; 
p' (AB) = 21 (00) p(AB) + 21 (10)p(A) p(B); (25) 
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(Aa) = 21 (00) p (Aa) + 21 (10) p(A)p (a). (26) 


To evaluate them we shall need, as in previous papers, the 


Lemma A). If Xn = 0X, + Ya (n=0,1,----) where lo) ai 
and if lim y, =y, then lim x, = y/(1—2) . 


noo no 


We now add the recurrence formulas for p’(A?) and p’(Aq) using 
(26). The result is that 


DE (Alia pe Uae 
and there is an analogous result for p™ (B). Thus 


pir) (A) == pen Ay, (Gi (a): 
p(B) =p (By; p™ (b) = BOO) 1 Ge ele): 


Next we apply Lemma A) to (26). Here a = 21 ® (00) = 21,(0) <1 
because of (19). Hence, considering that 21 (00) + 21 (10) 
= 1, it follows from Lemma A) that 


lim p) (Aa) — p) (A)p (a) ; 


and in the same way, since 2/,(0) < 1, it follows that 
lim p™ (Bb) =p (B)p© (by. 


n> 


Also, since 21 ©) (00) = 21,.(00) < 1 it follows that 
lim p™ (AB) = p(A)p(B). 


Next consider (24) and use now p™, p°"), and p = p: 


pint) (AaB) — 2](0) (4) (00) p™ (AaB) 
= 21 (01) p"(Aa) p(B) + 21 (10) - [p(A)p™ (aB) 
+ p(a)p™ (AB) ]. 


Asn ~ oo the right side converges toward 
p(A)p(a)p(B) (20 (10) + 210 (01) ] 
= [1 — 210 (00) ]p(A) p(a) p(B). 


Now 
1) (00) =1(00) +1(01) < 1,(0) <4. 


Hence the conditions of Lemma A) are again satisfied and 


lim p (AaB) = p(A)p(a)p(B) 


no 
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and an analogous conclusion holds for p’ (ABb), etc. Finally since 
21(00) < 1, we obtain the following result which covers the “general 
case” in which 1,(0) < 1/2,L(0) < 1/2. 

If (19) holds, then 


lim p™ (AaBb) = p(A)p(a)p(B)p(b). (27) 

2) Next assume 1,(0) = 1/2,1(0) < 1/2. From the definitions 
of 1,(0), 2.(0), and 1,.(00) together with equation (9) of I it follows 
that in this case 1(10) = 1(11) = 0, 1(00) + 1(20) < 1/2, and 
12(00) < 1/2. The recurrence relation (23) then becomes 


p') (AaBb) = 21(00) p™ (AaBb) 
+ 21(01) [p™ (Aab) p(B) + p™ (AaB) p(b)] (28) 
+ 21(02)p™ (Aa) p™ (Bb). 


Since 1,(0) = 1/2 it follows that p™ (Aa) = p® (Aa) = p(Aa) 

(n—0,1,----). In (28) 1(01) # 0 and 21(00) < 1. In the relation 

for p’“) (AaB) the coefficient of p™ (AaB) is 21(00) + 21(01) and 

this is, because of (9), less than one. Hence we have with p = p: 
lim p™ (AaBb) =p(Aa)p(B)p(b). 

We reason in the same way if 1,(0) < 1/2, 4(0) = 1/2. 
Finally if 1,(0) = 4(0) =1/2, and consequently 1,(1) =L(1) =0, 
it follows that 1(10) =1(01) =7(11) = 0 and 27(00) + 27(02) =1. 
If these last two values are both different from zero the recurrence 
formula for p°*» (AaBb) together with p (Aa) = p(Aa), p™ (Bb) 
= p(Bb), (n=0, 1, ----) shows that p™(AaBbb) > p(Aa)p (Bb). 
If 1(20) =1/2,1(00) = 0 the same limit is reached in the first filial 
generation, i.e. p™ (AaBb) = p(Aa)p(Bb) (n=0,1,----). If, final- 
ly, 21(00) = 1 we have complete linkage, and p™ (AaBb) = p (A@Bb) 
for n=0,1,----. Thence the theorem: 

If, for m=s=2 we have, in the “general case:”’ 


1,(0) =1(00) + 21(01) + 1(02) <4; 
1, (Oy = 1(00) + 21(10) + 1(20) <4; 


then, writing p = p we have for each allele 
p” (A) =p(A), +--+, p™ (B) = p(B), ++ (M0, 1,----)5 (29) 
lim p” (AaBb) = p(A)p (a) p(B) p(b) . (30) 


No 


If 1,(0) =1/2,1,(0) < 1/2, then 
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p™ (Aa) =p(Aa), p™ (B) = p(B), p™ (b) = pb), 
(n=—0 ,1,----); 
lim p” (AaBb) =p(Aa)p(B)p(b), (32) 


and similar results hold if 1,(0) < 1/2,1.(0) =1/2. 
If 1,(0) =1.(0) = 1/2, but 1(00) < 1/2, then 
p”™ (Aa) =p(Aa), p™ (Bb) =p(Bb), (m=0,1,----); (83) 
lim p™ (AaBb) = p(Aa)p(Bb). (34) 


Finally, if 1(00) = 1/2, all other |-values are equal to zero and, 
in this case of complete linkage: 


p” (AaBb) =p(AaBb), (n—0,1,----). (35) 


The corresponding limit results regarding the marginal distribu- 
tions p™ (AB), and p™ (Aa) of our distribution p™ (Aa@Bb) are the 
well known limit results of Haldane and Jennings. For Haldane’s s.d. 
(see section 3 of I) 1,(0) =1,(0) =1/6, hence the conditions (19) 
and (20) are satisfied and the limit result (30) thus holds for Hal- 
dane’s linkage theory. 


4, The Main Limit Theorem. (General Case.) 


We shall consider in this section the “general case” which corre- 
sponds to (30) and to the situation of the main theorem of Geiringer 
(1944). A complete theorem covering the various cases of “degenera- 
tion” of the s.d. (analogous to that in Geiringer, 1945) will then be 
established in the last section. 

For the sake of brevity we shall speak here of a “mixed” gamete 
(with respect to a certain locus) if with respect to this locus the gam- 
ete contains paternal as well as maternal heritage. The decisive condi- 
tion in the limit theorem of Geiringer (1948a) (m=1 , general “‘s”) was 
that 1(0) < 1/2, ie., that the probability 1-2 1(0) of mixed gametes 
be different from zero, while the decisive condition in the main theo- 
rem of Geiringer (1944) (s=1, general “m’”) was that all recom- 
bination probabilities (r.p.) c:y = 1 — 2li,(00) be different from zero. 
Let us write these conditions in the notations of the present paper 
where they apply to certain marginal distributions of our s.d. Re- 
membering (37) and (39) of I we first consider 


[(8-1).- (8-1) (g, oes ed bn et Be ge (e, pores Sm) : 


(e,=0,1), (G=1,----,m). (36) 
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This is an ordinary linkage distribution as in Geiringer (1944). From 
he...m(€1,**** a) We derive in the usual way by (m—2)-tuple summnia- 
tion the marginal distribution lj (e:e,). In our present notation this 
reads, fori=1,k=2, 


Lie (€1€2) == 1-1) (4-2) (8)... (8) (é, £2 0 «--- 0) = big Goh fe=)) (e, E2) } 


(ec, =0,1,2—0,1), oe 
and there are m (m—1) /2 such distributions. We also have 
UO) (2-00 (8) (e; ’ 0 ele: 0) == L (;); (ce, =0 ’ 1 ee 8), 
and there are m such distributions. 
It will be necessary that 
and 
lL: (0) S42 GS 1S em) (39) 


be satisfied. We shall see that these obviously necessary conditions 
are sufficient as well. Moreover we shall see (Lemma D) that (88) 
is a consequence of (39). It will be convenient to introduce a short 
name for 1-2/;(0). In analogy to 1-21;,(00) = ei, where cix is the r.p., 
put 1-27;(0) =e; and call c; “exchange probability.” Then (39) states 
that all exchange probabilities are greater than zero. 

Besides Lemma A) of the preceding section we shall need three 
more Lemmas: 


Lemma B). With the definition (39) of I we have forv=0,1, 
---- (Ss—2): 


Jr.) (¢,) =LM (e,) =F, (a=0,1,--- (sy) 


if, and only tf, 
[(0) (8)...-(8) = f, (8,) =% ’ (e,=0,1----8). 

In fact for m=1 this is the Lemma B) of Geiringer (1948). For 
general m it is no new theorem either, since ](9)(8)---(8) (g,) as well as 
J(2)---() (g,). are nothing but s.d.’s for m=1 and general “9,” - The 
equation which corresponds to equation (38) in Geiringer (1948) 
now reads 

i)(2,) = J, (8-2) (8) +++ (8) (e:) =4, rena US (40) 
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Next we need the: 
Lemma C). The inequality 
LOM (90----0) S 16979) (0.0 --++ 0) (41) 
holds if 
Cate] = pi ae (42) 


This is almost obvious, since all terms contained in /--- (00 ----0) 
are contained in 1) )--®) (0 0 ---- 0) as well, and if at least one of 
the inequalities (42) holds the second expression will contain terms 
not contained in the first. (Consider as an example l™® @ (000) com- 
pared with 1 ®@@ (000) which are given, symbolically by the expres- 
sions (J) + 1.) (lb + 24 + L)( +) and ( + 4) ( + 24 + b) 
(1, + 21, + L) respectively.) 

Let us return to the general recurrence formula (9) and the 
corresponding formulas for the marginal distributions which occur 
on the right-hand side of (9). Consider first all marginal distribu- 
tions of order two, such as (avoiding double subscripts) p™ (A?), 
p™ (Aa), ----p™ (AB), etc. We have with 


[ (8-2) (8)... (8) (e) —y He) (s) : 
pn (A?) = 21,0) (0)p™ (A?) + 2L° (1) [p(A)]?. 


Here, because of (39) and Lemma B) the coefficient of p™ (A?) is 
less than one, hence, because of Lemma A), p (A?) > p(A)?. We 
also have 
po) (AB) + Bho) (00) p> (AB) 
Ge) (10) D(A) OCs): 


Here 21,2.) (00) < 1 because of (38). Using Lemma A) we see 
that p™ (AB) > p(A)p(B). 

Next we consider the recurrence formulas for the marginal 
distributions of order three, such as p) (A?) , ,p™ (A2B) ,--- 
p™ (ABC) ,---. For p™ (AS) the convergence ret [yp (A) }? fol- 
lows because of (39). We also have 


pt) (AB) = 1, (2-2) ) (00) p™ (A2B) + .--++. 


Now 1,2.) (00) < 1.6 (00) because of Lemma C), and this 
last value is less than 1/2 according to (38). Hence because of Lem- 
ma A), and of the convergence of the distributions of order two we 
see that p™ (A?B) > »(A)*9(B). Also in 


p21) (ABC) = Qljog (8-2 (8-2) (8-2) (000) p™ (ABC) +..-. 
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we have 1,238) (1) (8-1) (000) < byes * 048) (000) = 1.1) (*-) (00) and 
since this last one is less than 1/2, the coefficient of p™ (ABC) is 
less than one and by means of Lemma A), it follows that p\” (ABC) 
> p(A)p(B)p(C). 

We continue in the same way. The coefficient of p™ (A*) in 
the recurrence formula for p(» (A) is 21,(-® (0) and this is less than 
one because of Lemma B) and (39), and the coefficient of a term like 
p (A*B8CY) in the respective formula is 


Qliog 8-® (8-B) (8-7) (8)... (8) (000) < 21. (8-1) (s-) (00) ah ‘ 


because of Lemma C) and (38). Hence we may each time apply Lem- 
ma A) and the convergence toward the product of all marginal dis- 
tributions of first order follows. 

Before formulating a final result we prove: 


Lemma D). If a recombination probability ¢,; = 1-21,;(00) van- 
ishes, then the two corresponding exchange probabilities c; = 1-21, (0), 
and c; = 1-21,(0) are zero as well. 

First let us prove this for m=2 and general “‘s.’”’ We have, sym- 
bolically, 


262(00))— 2[, (es —1) hh La? 
Now if 
21,2(00) =1, 
and since the sum of all /-values equals one, 
1 [hi Sh ee 12, 


we see, considering the symmetry relations for the s.d., that the co- 
efficient of each single term to the right in the second of these equiva- 
lence relations is greater than of the corresponding term in the first 
one with the only exception of 1(00). Hence, if we subtract the smaller 
sum from the larger one it follows that each /-value equals zero ex- — 
cept 1(00). But if 21(00) = 1, ce = 0, then 4.(0) = L(0) = 1/25 
and this proves our statement. 

Now consider general m and s. Corresponding to the ith and jth 
factor we form the marginal distribution of order two for general s. 
For i=1 , j=2 this gives 

J(0) (0) (8)---(8) (gg, & O-++- 0) = Lin © (€, &2),, (a =0, 1,-::-s). 


This distribution has all the properties of a s.d. for m=2 , and gen- 
eral s. Hence for this distribution Lemma D) holds, which means if 
Lo?) (00) = 1.(00) = 4, then 1,(0) = E(0) ==" 1/2 Here 
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lio) 6) (ge) = le (ee) has been derived from hh» (e:¢2) by an 
2.- (s—1)-tuple summation, and J, (¢,) from the same distribution by 
an s-tuple summation. But these same distributions are also the mar- 
ginal distributions of our original s.d. with arbitrary m . In fact 


bia (é; @) = [eevee 0 ----0), (250, 1,20, 1), 
L, (€1) = [(0)(8).:..(8) (9, ,O---: Oh. (e,—0 seal kes +202 $) . 
bP (a5) HT). (0, e,,--+-0), (e.=0,1,---- 8). 


Hence the Lemma is proved for the subscripts 1 and 2 and therefore 
for any pair of subscripts. 

Let us note that conditions (39) are not only sufficient but neces- 
sary as well, considering the limit result we have in mind. This will be 
clarified in the following section. However, using this fact in the for- 
mulation of our theorem we obtain the following general result: 


Theroem I. If, and only if, conditions (39) hold, or in words, if 
and only if, the m exchange probabilities ec, = 1-21,(0) are positive 
the distribution of gametes converges toward the product of thems 
marginal distributions of the single genes in the initial generation. 
This holds under random mating and chromosome segregation. 

From this theorem follows a result for the d.zy. by means of 
equation (4) of I. 


5. The General Limit Theorem. 


We still have to consider the various possibilities connected with 
a “degenerated” s.d. The “general case” studied in the preceding 
section gave the result: 


pi) (a,, vies Uys tts Amy +++ + Ams) > p(a,,) -+++D (Ans), (43) 


if, and only if, no exchange probability (and consequently no recom- 
bination probability) is zero. If, on the other hand, all r.p.’s (and 
consequently all exchange probabilities) are zero then: 


DP (Air °+** Gms) =P (dar +++ Ome), (N=0,1,----). (44) 


In these last formulas we again used the notation p® =p. Between 
these two extreme cases are the various cases of groups of linked 
genes. 

Let us first consider an example which will exhibit the various 
possibilities. Assume m=4 , s=3 , ¢, =c.=0,¢, #0, ¢, #0, and all 
six r.p.’s different from zero. Avoiding double subscripts we first con- 
sider the marginal distribution of the three genes at the first locus. 
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The recurrence formula then reads 

DP") (GA2d3) = 21, (0) (a,do43) + +--+, 
Since c, = 0, ie., 21,(0) =1, it follows that 

D™ (AyA2A3) =P (Ayo), (N=0,1,----), 


which means these three alleles are completely linked. In the same 
way it is seen that, because of c, = 0: 


PD OCG020,) =p (0.0.05), -(n=0'; 1, >="). 
On the other hand consider the recurrence formula for 
p'™*) (A,0,03b,b,b3). 


The coefficients at the right-hand side of this recurrence formula are 


the values of 1,5) (e, 2), (e; =0,-+-- 3). Let us write for the mo- 
ment /(e, e2) instead of 1,2. (e, e2). There are 16 such terms. Since 
L(1) = 0, 4(1) = 0, and the corresponding symmetry relations 


hold, it follows that of these 16 terms 12 are zero and 1(03) = 1(30), 
1(00) = 1(83) are the only ones not necessarily zero. Next we see 
that 1(03) =1(80) # 0 because, as can be easily seen, ¢,. would other- 
wise vanish. If 1(00) + 0 we find, as usual, by means of Lemma A) 
that 

Pp” (A,4,0;b:b2b3) > P(A 1Aot3) P (b1b2b;), 


and, if 1(00) = 0, it follows that for all n: 
p (A,d2030,b2b3) = p (a1d.a;) p (b1b2b3). 
On the other hand, since c; # 0, we see easily that 
Dp” (€,6,63) > D(C) DP (C2) p (Cs), 


and similarly for the fourth locus. 
If, finally, we consider the joint distribution of all four loci it 
may be seen by analogous considerations that 
lim p\” (@,4.0; «+++ d,d2d;) =p (a,a.0;)p(b,b2b;)p (C1) «+++ p (ds) ° 


No 


As a second example consider the same situation but with ¢c.=0. 
It then follows for the linkage distribution I(« 62) =e" (e: €2) 
that 1(00) = 1/2, and consequently 


pi) (@,a,a;b,b,b;) = 1b (@,4,0;b,b,b;3) , (n=0,1,--- ‘), 
and for the joint distribution 
lim p\” (a, +--+ ds) =p (a, +--+ bs) p (G1) -++- p (ds). 


nro 
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If, next, not only ¢,. but also ¢,; is equal to zero it follows from 
the “triangular relation” for the r.p.’s (see Geiringer, 1944, page 34) 
that ¢,; = 0, and from Lemma D) that c, = ¢ =c¢c, =0. Assume 
c, > 0, then necessarily C144 > 0, Gos > 0, Cag > O, and we obtain 


lim p (a, --+- ds) =p(di--++ ¢s)p (di) p (de) p (ds) . 
In fact, according to the above results we may say: 

If ec, = 0 there is complete linkage between the s genes which 
correspond to the ith locus. If c,; = 0, and consequently c, =c; = 0, 
there is complete linkage between the 2s genes which correspond to 
these two loci. 

If we want to completely describe the various possibilities we have 
to introduce the concept of maximal groups of completely linked fac- 
tors, used first in Geiringer (1945). Consider the m numbers 1, 2, 

--m. There are t = 0 disjoint groups G,, G2, -+:- G; each contain- 
ing uw; 2 1 of the m numbers 1,2, ---- m (@=1, ---- t) such that 
within each group all r.p.’s are zero but between any two such 
groups at least one r.p. is-different from zero; this holds if u, 2 2. 
In addition we speak of a “group of size one” (uw;=1) if the corre- 
sponding ¢;—=0 (the exchange probabilities which belong to numbers 
contained in a group of size > 1, are a fortiori zero). The remain- 
ing numbers which do not belong to any “group” and for which there- 
fore the corresponding c; > 0, we call “free elements.” We may de- 
scribe the situation as follows: Within each group no recombination 
takes place, but there is at least one r.p. different from zero between 
any two Groups G; and G,, and between any “free” element and any 
other number (whether “free” or in a “group’’). If, for example, for 
m=i we know that¢, 0 ¢:. 4 0, Gis = Ca — C4. — 0 ter 
angular relation ¢.; = ¢ + C.; shows that ¢., — 0: hence (4, 6, 7) 
form a group, but this is not a “maximal” group. Moreover (1, 2) 
form a group, since ¢:. = 0; but in addition c,, = 0; thence, because 
Of Cos = Cro + Cig, Cog = 0, and continuing that way we see that 
Cra = Crp = Chz = Coq = Cog = Cor = O and the five elements (1,2,4,6,7) 
form a maximal group of completely linked factors while “3” and 
“5” are “free elements.” 

Now assume that G,, ---- G; are the disjoint maximal groups of 
completely linked factors. Then for the gametic probability corre- 


responding to each such group (marginal distribution of order s - 2;) 
we have 


Dat” = De, (t=1,----t), (n=0,1,----), (45) 


which means the complex of s - 7; genes corresponding to G; forms a 


a 
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group of complete linkage and is conserved through the generations. 
This follows in the usual way by writing down the recurrence for- 
mula for Pa." where on the right-hand side the coefficient of De” 


equals one. Next consider the joint (marginal) distribution of two 
groups G, and G, , namely Deas We see in the same way as in the 
examples, using the fact that no r.p. between the numbers of G, and 


the numbers of G, equals zero, that the recurrence formula is of the 
form 


pee eg. De a, + BPM a. (Os 
with a + 6 =1and £ > 0, and consequently either 


DOs @: = De, Des (Oe Lato) SOF Pra a = Da Da, : 


Hence certainly lim Pg a, De, Pc,. Following up these ideas we 
no 


see that 


lim DG Gaal, = De, HAI tet De, : 
No 


On the other hand the s genes which correspond to a free element are 
not linked since the corresponding exchange probability is different 
from zero. If 7 is the number of such an element we have 

lim p™ (a, aerete a, ) =p(a.) op (G, Mee 


No 


We may formulate these facts in the following 


Theorem II. For a 2s-ploid of r alleles and m loci consider the * 
heredity assumptions of this paper (random mating, chromosome seg- 
regation, etc.). Denote by Gi,----G, (9St= m) the peace groups 
of completely linked genes,” each group containing ig eae numbers 
(1 ---- t),. (a te + U1 =m), end by), ead k the remain- 
ing f free elements (0 = f S m—u) with non-vanishing exchange 
probabilities. Then: 

1) The complex of (s- ui) genes corresponding to each group 


G, remains permanently linked, while with respect to the free ele- 
ments the distribution of each single gene only is preserved 


n —- ps shaveney ae 46 
PNG Sie yan pe at Dat (n=—0, 1, as (46) 
9 (jc) =P (Ajo), +77 1 (Ao) =p (Uo) ; 


(47) 
(n=0,1-:*), (o=1,----8). 
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2) The joint distribution of allm-s genes converges toward a 
limit such that the original marginal probabilities corresponding to 
each maximal group as well as the f - s original marginal probabil- 
ities of the remaining single genes are independently distributed. 

If f=m , which means if there are free elements only, our theo- 
rem reduces to that of the preceding section. For s=1, it gives the 
general result of Geiringer (1945). (If s=1 the new concept of a 
“oroup” of size one with vanishing exchange probability has no re- 
alization.) If m=1 and general s the exchange probability, c, with 
respect to the one locus is either zero or not and accordingly (41) or 
(43) of Geiringer (1948a) holds. 

As an example consider m=9 , s=4, and assume we know that 
Ce = 63 Ce = 0 6 6a = C505 0 DUT CeO Cre eee 
It follows that G, = (1,2,3,4), G. = (5,6,7), G; = (8) are the three 
maximal groups and “9” is a free element. We then have: 


lim p" (Qyq + +++ yg) == PD (Ayy +++ Ugg) DP (M51 +--+ Mra) 


noo 


‘DP (Qgi+**+ Asa) DP (M1) +++ DP (Ana). 


In formulating our main result we may as well start with the 
consideration of the free elements and determine the ¢ disjoint maxi- 
mal groups of complete linkage for the remaining (m—f) numbers. 
This gives: 


Theorem III, Denote by c,, @ , -++- ¢¢ the non vanishing exchange 
probabilities (0 =f S m), and by G,, Ge, ---- G; the maximal groups 
of completely linked genes, each containing u, 2 1 of the remaining 
(m—f) numbers (u, +---- u, = m—f). Then under the hypotheses of 
this paper: a) The original (marginal) probability of each of the 
(m-s) single genes is preserved through the generations. b) Within 
each group the original (marginal) probability of the whole group is 
preserved. c) The joint gametic distribution (m - s-variate distri- 
bution) in successive generations approaches a limit where genes and. 
groups of genes are independently distributed. 

From these theorems we derive corresponding limit results for 
the distributions of zygotes by means of equation (4) of I. 
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A geometrically constrained enzyme system is presented which is 
capable of driving a substance, which acts as a coenzyme to one of the 
enzymes of the system, against a concentration gradient. The energy for 
the process would be derived from the breakdown of the substrate for 
the enzyme-coenzyme system. If auxin-like substances act in the man- 
ner of a coenzyme in plants, then the mechanism presented might serve 
as a model to explain the polar transport of such substances in ‘plants. 


In the following paper a mechanism is outlined which could trans- 
port polarly, against a concentration gradient, a substance that acts 
as a coenzyme to some enzyme constrained in the transporting cells. 
The energy for the work done against the concentration gradient 
would be derived from the breakdown of the substrate of the enzyme- 
coenzyme system. 

If such growth stimulating substances as indole-3-acetic acid and 
the so-called “auxins” in plants act as coenzymes to some enzyme 
fixed in the plant cells (Commoner and Thimann, 1941; Sweeney and 
Thimann, 1942; Bose, Dutt and Thakurta, 1944-1946), then the mech- 
anism to be described might serve as a theoretical model to explain 
the polar transport of these substances in plant coleoptiles. Although 
this theoretical model was conceived with reference to the problem 
of polar transport of “growth hormones” in plant coleoptiles, it should 
be pointed out that its possible application is by no means restricted 
to this particular case. 

This type of mechanism bears a similarity to a more general 
one, which has been discussed previously by J. Franck and J. E. Mayer 
(1937) without any reference to this particular problem, and also 
to the mechanism which has been proposed for active uptake of glu- 
cose into the gut by means of phosphorylation and dephosphorylation 

Hober, 1945). as 
isi the bade theory, the coenzyme-like substance, which m 
brevity shall simply be designated as “J,” is transported the iepet 
of each cell by protoplasmic streaming, but transported polarly acy 
the boundary regions from one cell to the next by a “pumping” mech- 
anism to be discussed presently. The protoplasm at the end regions 
is, in this model, to possess the mechanism for polar transport. 
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If the model is to be applied to the specific problem of auxin 
transport, then the transporting cells, indicated schematically in Fig- 
ure 1, would be a series of parenchymal cells of the coleoptile. 


DIRECTION OF TRANSPORT ————> 
re 
ae) Cert poo 
iD) | ee) | (CEES KE [GES 
’,. / é 


PUMPING PUMPING PUMPING 


MECHANISM MECHANISM MECHANISM 
(Bounpary REGION) (BOUNDARY REGION) (SOUNDARY REGION) 
FIGURE 1 


The streaming protoplasm acts in the manner of a “conveyer 
belt” and the pumping mechanism of the end regions serves to un- 
load I from the conveyer of one cell and load it onto the conveyor 
of the next cell. 


4-8) 
FIGURE 2 

Figure 2 shows a highly magnified diagram of the boundary re- 
gion, whose relation to the rest of the cell is indicated by Figure 1. 
It is this mechanism, whose geometry is illustrated in Figure 2 (a 
discussion of the kinetics will occupy the major portion of the remain- 
ing part of the paper), which is conceived as driving J against the con- 
centration gradient. 

Let us consider the mechanism of Figure 2. Suppose the proto- 
plasmic layer comprising this boundary region is divided into three 
layers—2, 3, and 4. Let the entire system be bounded on both sides 
by membranes (1-2) between region 1 and the layer 2 and (4-5) be- 
tween layer 4 and region 5. Layer 2 is distinguished from the other 
layers of the system by the fact that it contains imbedded in it an 
enzyme, E', for which I is the coenzyme. The protoplasm between 
the two membranes which bound the system is visualized as being 
gel-like in its structure so that the enzyme, EF’, is geometrically con- 
strained to layer 2. Layer 3 is neutral insofar as the particular set of 
reactions under consideration is concerned and enters into the present 
problem solely as an impedance to diffusing substances. Although 
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layer 3 is treated as merely a homogeneous gel in this discussion, in 
the plant it would probably consist of a relatively impermeable lignin 
and cellulose wall pierced here and there by threads of protoplasmic 
gel through which lighter molecular weight substances could diffuse 
from one cell to the next. As long as layer 3 is not completely imper- 
meable, however, its precise internal structure is of no concern to 
the present problem. Layer 4 contains, imbedded in its gel, an en- 
zyme A. 

One of the simplest and most common types of general reactions 
involving an enzyme E and a coenzyme J with some substrate S is: 


ky 
EFI EI: 
k, 


Ks 
EI+S=[ EIS; (1) 
ks 


ks 
EIS > E' + IB + Other products. 


The substance JB then reacts with A , the complex breaking down 
to release J and A again: 


Al LB ALB 
AIB SASSI EB . 


In our model, reactions (1) take place in layer 2, IB diffuses 
from layer 2 to layer 4 where reactions (2) occur, releasing J. Since 
I is consumed in layer 2, the concentration of J there will normally 
be less than in 1 and also less than in that part of region 3 immedi- 
ately adjacent to 2. Thus J will diffuse into layer 2 from both 1 and 38. 
The quantity J is being released in layer 4, however, so that the con- 
centration of J in 4 will be higher than in the immediately adjacent 
regions of both 3 and 5, provided certain conditions, which will be 
discussed later, are met. From the preceding considerations of a 
qualitative nature it is apparent that part of J will diffuse back across 
layer 3 to 2 and part of it will diffuse into 5. 

When the system is operating in a steady state, J will be enter- 
ing from 1 as fast as it is diffusing into 5. The net overall effect 1S 
that I is being transported from 1 to 5 and if the concentration of J 
in 5 is greater than in 1, this will amount to active transport, i.e., J 
is being transported against a concentration difference. The assump- 
tion has been made that the bounding membranes are permeable to 
I but: impermeable to the larger molecules of JB, so that IB can only 


(2) 
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diffuse through layer 3 from 2, where it is formed, to 4, where it is 
consumed. It is not absolutely necessary to the working of the model 
that the membranes be completely impermeable to JB, but it simpli- 
fies the mathematical analysis to assume complete impermeability to 
IB. 

The energy for the pumping action comes from the overall reac- 
tion, which is the breakdown of S into “other products” and B. 

In the subsequent discussion, we shall introduce the following 
notations: 


I,= concentration of J in 1; 
I, = average concentration of J in 2; 
I,= average concentration of I in 4; 
I; = concentration of I in 5; 
62 = thickness of layer 2; 
6; = thickness of layer 3; 
6, — thickness of layer 4; 
IB, = average concentration of JB in 2; 
IB,= average concentration of JB in 4; 
V2 —= average velocity of reaction (1) in layer 2; 
v,— average velocity of reaction (2) in layer 4; 
A=) — le 


Care must be exercised in making the kinetics formulations on 
the assumption of steady state conditions, because we have said that 
the last stage of reactions (1) and (2) proceeds irreversibly. If J. 
is made sufficiently large so that 


b2V5 a OsV4 max 9 (3) 


then a steady state is impossible, since the amount of JB in the sys- 
tem will continue to increase without limit. To insure against this 
defect in the formulation, one of two alternatives is available. First, 
the last stage of reaction (1) could be made reversible, but then some 
assumption would have to be made about the nature of the products. 
Secondly, one could assume conditions such that v. was the limiting 
reaction, i.e., such that 


O2V2 max << O4V4 max ° (4) 


The latter assumption will be chosen because it is less drastic 
and less restricting than the assumption of some specific nature of 
the products. More will be said of this after the kinetics formula- 
tions have been made. 


pi Ne 


ay 
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A flow scheme for the various substances involved in the mech- 


FIGURE 3 


anism is indicated in Figure 3. The broken arrows labeled v, and v, 
indicate chemical reactions while the full arrows indicate directions 
of diffusion. Having now given a qualitative description of the theo- 
retical model, let us proceed to a quantitative investigation of its 
properties. 


The Kinetics Formulation. 


Let k,, k., k., k» and k; be the reaction velocity constants as 
indicated for reactions (1). Then we have for steady state conditions: 


V2 = hy (B) 1) — ka (ED) ; 
Oe (hi) (yt he (BIS) — h(E) ke (Ee]) (S)’; (5) 
0 = ka (ET) (S) — k2(BIS) — k, (BIS). 


If EH, is the total amount of enzyme per unit volume in iar 2, 
then 


Oa ie lS (6) 


Solution of equations (5) and (6) gives 
(E.) U2) (S) kikcks 


= eG) 
02 Teales (In) (S) + Tey (Vig + Ba) (Ls) + Bakes (S) + fea (a + hs) 
Similarly for the reaction (2) occurring in layer 4, we have 
(A,) UB.) ks (IB,) W 
y= = ey (8) 
Wa ates (Bi) ae 
GRE Se ( =) 
where 
WwW = Ks (A,) ’ 


Kh heh 


Let us now see what assumption (4) involves. For a fixed value 
of S 
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Ey k.k3S 
Uses 
We keSurig eo. (Kane Ba) 


as I, > large, (9) 


or, if both J, and S become large, 
V_ > E’yks = Ve max (10) 
Then the condition under which the steady state formulation ap- 
plies is that 


Os 
Ek; << Fe KAS: 


FIGURE 4 


Figure 4 illustrates graphically v. as a function of J, for vari- 
ous concentrations of S (full lines). On the same set of axes and to 
the same scale is also shown 


as a function of JB, (broken line curve). 
From assumption (4) it follows that H,k; must always be less 
than 


and so, since for steady state 


64 
V5 << Se Vs; 


the system will usually operate in the region where 
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is approximately linear. This approximation does not enter into or 
affect the final result since v, completely fixes v,, but it does some- 
what simplify the algebra involved in the elimination of v4, from our 
equations. 


Diffusion Relations and Boundary Conditions for Steady State. 


Let D,= diffusion coefficient of J through 3 ;: 
D, = diffusion coefficient of IB through 3; 
V = flux of J through the system (amount of transport), 
expressed in gm cm” sec; 


h, = permeability of membrane 1-2 to I; 
h,—= permeability of membrane 4-5 tol. 


From condition of steady state, the rate of consumption of I by 

reaction (1) must be the rate at which J diffuses into layer 2, so that 
i, Sa Ib 

620, = h, (1, — Ie) + D, 5 ° (11) 


3 


Similarly, the rate at which J is liberated by reaction (2) in layer 4 
must be equal to the rate of diffusion of J out of 4, so that 
I,—I 
6.0, =h AL, —Is) + Di 7 =, (12) 


3 


Since the membrane 1-2 is impermeable to IB, the rate of produc- 
tion of JB by reaction (1) is equal to the rate of diffusion of IB across 
layer 3 and the rate of consumption of JB there by reaction (2), hence 


Lo 1D. 
02V2 — D, a (13) 
63 
—IB 
O4V4 — 10); es 5 i. \e (14) 
3 


If the total amount of the various forms of J in the system is to re- 
main constant, then the inflow of J across one of the membranes must 
be equal to the outflow of J across the other membrane, so that 


V=h(4—1). (15) 


From the set of equations (11)—(15) we obtain, after some alge- 
braic manipulation, the expression 
fee D, D,h, oi D,hy of O,Noh, ; (16) 
* 8085 5, Sshah, 
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The form of (16) may be somewhat simplified if we put 


st hat he ; (17) 
hihe 
so that (16) takes the form 
iH ae ro) 
pes Dud i VAD Hei eee (18) 
6» 63 62 03 
The Solution for V. 


If we solve (15) for I, explicitly and substitute this into the 
kinetics formulation (7), putting 


we find 


Vo = B. (19) 


For steady state operation, the two expressions (18) and (19) for 
V2 must be satisfied simultaneously. 

For convenience in manipulation and to make more apparent the 
form of the relations, put 


M= mae (20) 
62 03 
Nets (21) 
62 0s 
Q=(S+K); (22) 
a KS eK ey 
ky 
R=KES. (24) 


The behavior of V can best be examined graphically (Figure 5). 
If v, is plotted against V, equation (19) gives rise to a two branched 
curve, with a different pair of curves for each value of I,. The ver- 
tical asymptote is shifted to the right along the positive direction of 
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FIGURE 5 


the ve axis for increasing values of I, , while the horizontal asymptote 
remains the same for all values of I,. The vertical asymptote of the 
pair of § curves will obviously always lie in the right half of the co- 
ordinate plane, since J, > 0 and Ph./Q > 0. 

Equation (18) represents a family of straight lines, a, of con- 
stant slope but with different intercepts on the v. axis. Each value 
of A gives rise to a different « line. The intersections O and O’ of the 
a and § give the two values of V that arise from simultaneous solu- 
tion of (18) and (19). It will be shown, however, that only the inter- 
section O has any physical meaning. It will also be shown later that 
this intersection O represents a stable configuration point for the 
system. 

It should be noted that the straight line y represents a bound for 
all the a lines. That is to say, all the a lines must lie above y. The 
line y has the same slope as the a lines and intersects the V axis at 
I,h.. The reason for arises from the following physical considera- 
tions. 
For an a line to cross the V axis at a value larger than [hz , it 
would be necessary for it to meet the condition that 
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{ 
ee STI (25) 
N 
or that 
D,A 
Gs SSF PSS SS > 1G , (26) 
DH 03 


which is tantamount to saying that in the pure diffusion case, i.e., 
with no driving reaction (see page 234 of this paper), the diffusion 
across the entire system is greater than the flow across the bounding 
membrane (1-2). This is obviously wmpossible. 

Further information can be obtained from Figure 5. The limit- 
ing value of A against which the system can pump, for infinitely 
large I, is 

R82 63 ks HS 


Aa (27) 
MQ D,(S + K) 


The limiting or maximum transport, for infinitely large J, , corre- 
sponding to any value of A is given by 


(= am )— V h nied | Bette 
ae eae —— V max, Where — — Ss 5 
Q N N 


or 


ees Tee ds0s ea 
(S+:K) ~ 6,8, 1 (D,H ¥8;) _ 


Analysis of the Non-Steady State. 

Up to now in our investigation of the theoretical model, which 
has been postulated, the steady state has been assumed in all of the 
formulations. Two roots in the solution for V were obtained, giving 
two values of V corresponding to a given J, and 4. We shall now 
study these configuration points with regard to their stability, to as- 
certain whether they represent stable dynamic equilibria and which 
root would represent the actual steady state configuration point of 
the physical system. 

For layer 2 the rate of change in the amount of J in the layer 
is given by the expression 


Vmax: (28) 


dl. dD, 
Og ase te) = 7% (I, — I.) == Oa0s, (29) 


and for layer 4 by 
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Cay 4 a 5 4 . 


Since all of the layers will be thin, the approximation will be 
made that the gradient of JB across 3 always adjusts rapidly enough 
so that we can consider the flow of JB into 3 to be equal to the flow 
of JB out of 3. This approximation has also been made for I, in 
writing expressions (29) and (80). Since there is practically no 
delay in the adjustment of v. to new concentrations of IB, we can 
write as before 


O2V2 = b4V4. (31) 


Introducing equations (7), (20), (21), (22), (28), (24), and 
(31) we can now write (29) and (30) in the form 


i Bice OF R De 
Su fe ae Pe ere 
dt ass Oo 03 QI. =F EP. Oo 03 
dl, pmo ley a. ) fe a a) 
—*——* “43 zh: ins eT eeere(88) 
bash = (s+ A) CS ray EE et oe seas 


Then the curve dl./dt = 0 will be given by 


ha tO) ee 
(he b3+ o>): 


R_ he daly 
D, 
FIGURE 6 
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9 O28 ho 0 ae ID), Os 6,R fe 
ahs 1+ ( — =) 1. +( ) =I,, (34) 
D D, (QI, + P) 


1 BE 


and dI,/dt = 0 will be given by 


( D, i -( 5, dR Tee lige: oe) 


Fee rf ee , = I,. (35) 
h, 6; + D, eds Ds (QI, + P) (A, 63 + D,) 


If J, is plotted as a function of J, as in Figure 6, it may be seen that 
the curves representing (34) and (35) have only one point of inter- 
section for positive J, and I,. That this is true may be verified by 
noting that the slope of both curves approaches a constant value with 
increasing J, and the slope of (34) is greater than (35) for all posi- 
tive I,. Thus there is only one steady state root corresponding to 
dI./dt = 0 and dlI,/dt = 0 that is physically significant, i.e., not in- 
volving negative concentrations, that arises from the solution of (34) 
and (35). It may be seen from these equations that for J, greater 
and I, less than the values lying on the two curves, dI,/dt is less than 
0 and dl./dt is greater than 0. This would mean that any point in 
that region would move downward and to the right with time. With 
each point it is possible to associate two vector components, one along 
the direction of the J, axis and one parallel to the J, axis. A similar 
procedure may be carried out by examining the values of the two 
time derivatives for points of each of the other regions bounded by 
the two curves. This is done in Figure 6 and the longer arrows indi- 
cate the resultant direction of motion of points of the various regions 
of the I, vs. I, plane. Figure 6 shows that this point given by the 
intersection of the two curves is a stable point, since regardless of 
what the initial values of J, and J, are they will, with increasing time, 
approach those values corresponding to the point P. 

Solving (34) and (85) analytically to eliminate 7, and obtain 
a solution for J, in terms of J, and A , we have 


Ped “ Var ee AGP. 


where 
re Oo 032 Rh, 37) 
(hahsd:? + hi63D, + hyd3D,)” ( 
and 


(hoh.63? + hib;D, + h.d;D,) 
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The constant C is always positive, since when A is negative |A4| < J,. 
All other terms in C are positive. The quantity b will always be posi- 
tive, since it involves only positive quantities. Similarly we have seen 
previously that P and Q are both positive. Thus 


JS c) += o(=—* c) (39) 
Q Q Q 

so that only the root obtained by choosing the positive sign of the 

radical will have any physical significance, i.e., 


Paso PEO 2A GP, 
mE 


Examination of the case when A = 0 and equation (38) reduces to 
C = 1, shows that from (15) and (40) we have 


V=h,(1,—1.) ° 
(41) 
IP a= ip (0) (P20)? 
= = I asees ee We ae —-{/, + ——— ], 
ha( I, ely + 20 :,/ Q Q? ) 


and as I, > large 
Dae (P — b) 
Vie [ n-an ean (1+ 45> +) 


ie hb 1R 
Do met ti Qe 
which we see corresponds to the intersection O of Figure 5 for these 
conditions. 
If, on the other hand, we choose the root giving a negative value 
of I, , we obtain by making 4 = 0 and J, large: 


V=hl ga ele 
Laas 241 Q 5) 


which, inspection of Figure 5 will show, corresponds to the operating 
point O’ under these circumstances. Hence we see that O has no 
physical significance and the point O of Figure 5 is the operating 
point of the system. 
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The Pure Diffusion Case. 


Since reference has been made earlier in this paper to the de- 
generate case when there is no reaction, but merely diffusion across 
the system, it shall be derived at this point for purposes of reference. 
We can write immediately that 


he (I, — Ie) =h,(1,—Is) > (42) 
“ (I, —L,) =h,U. — 15) ; (43) 
V=he(h — Ia); (44) 
from which is obtained upon elimination of J, and J, that 
—D, A 
ee eee (45) 
D,H '+ 6; 
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BOOK REVIEW 


JOSEPH SCHILLINGER. The Mathematical Basis of the Arts. 1948. x.+ 696 pp. 
New York: Philosophical Library. $12.00. 


A mathematical theory of art is of interest both to the mathematician and 
i the mathematical biologist. In his pioneering work the late George, Birkhoff 
(1983) made an attempt to evaluate the effect of different types of, symmetry, 
order and other quantitative elements of a visual or acoustical pattern upon . 
its pleasantness. His approach was a purely formal one. He postulated, in a 
somewhat intuitive way, that certain geometric or temporal characteristics of 
a pattern possess a certain arbitrary degree of pleasantness or unpleasantness. 
Taking then a number of concrete patterns he evaluated, on the basis of his 
postulates, their relative aesthetic values. The importance of Birkhoff’s work 
is not diminished by the fact that experimental determinations by R: C. Davis 
(1936), who used the rank order method, of the actual relative aesthetic values 
of some of Birkhoff’s patterns decidedly disagreed with Birkhoff’s predictions. 
These experiments show that some of Birkhoff’s postulates are inadequate. The 
fundamental method of approach remains of interest, and the very fact that R. 
C. Davis’s experimental work was suggested by Birkhoff’s theoretical studies 
gives the latter at least the value of a good working hypothesis. 

A different approach is made by the mathematical biologist. He begins with 
a set of assumptions about the mechanism of the central nervous system, and 
of what processes in it correspond to a feeling of pleasantness or unpleasant- 
ness. On that basis he then develops a theory of perception, in. particular of 
visual patterns. The assumed nervous mechanism indicates what geometrical 
features of a pattern are pleasant and to what extent. The aesthetic value of a 
perceived object is thus considered a function not of the object only, but of the 
perceiving individual. Calculation of relative aesthetic values of some of Birk- 
_ hoff’s polygons, which were experimentally studied by R. C. Davis, led to quite 
a surprising agreement between theory and experiment. Subsequent- theoretical 
and experimental work further improved this agreement (Rashevsky, 1948). 

In view of the above, it is natural that an attractive looking volume of almost 
700 pages entitled, “The Mathematical Basis of the Arts,” should catch the eye 
and excite the curiosity of the mathematician and of the mathematical biologist. 
It seems, however, fairly safe to say that both the mathematician and,the mathe- 
matical biologist will be sorely disappointed in this book. This in itself, the re- 
viewer hastens to remark, does not imply any reflection on the quality of any 
book. 

Joseph Schillinger’s book undertakes a detailed classification* of different 
numerical relations which occur in various arts. Possible numerical sequences, 
ratios, symmetries, etc., are painstakingly classified. However, the reader will 
not find any mathematical argument, proof or deduction, as usually understood 
by the pure or applied mathematician. No discussion is given of the possible 
quantitative connection of the different numerical] relations with the -aesthetic 


*Italies by the reviewer. 
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value of a work of art. In no sense does the book present a mathematical theory 
which could be tested experimentally and which would explain anything. 

Such a classificatory study may be very valuable, perhaps with the view of 
being used as a sort of empirical material for later theoretical studies. It would 
have been, however, more advisable to call the book, “Numerical Basis of the 
Arts,” or, still better, “Numerical Relations in the Arts.” The prevailing lay- 
man’s opinion notwithstanding, numerical and mathematical are not synonymous! 

We shall illustrate the above statements only by two examples, which are 
rather typical of the whole book. 

On page 64 the author begins an enumeration of different “Art Forms.” 
Thus Art Form Number One has two aspects: (a) art of audible sound* (mu- 
sic); and (b) art of audible word* (poetry). Art Form Number Two is the 
“art of touchable mass.”* Number Four is “smellable odor.” In this way eighteen 
“Forms” are established. 

Then on page 75 we read that, “A complex homogeneous form may consist 


of combinations by two or three elements .... ,” and the possible combinations 
are tabulated. The tabulation takes 3 pages, of which the following is a sample: 
“By TWO 


5+8;5+9;5+10;5+4+11;5+12;5+13;5+14;5-+ 15; 
5 +16;5-+17;5-+ 18. 


S RH 

We see 3 pages of numerals, but hardly anyone would call this mathematics! 
Another example may be found on page 344, where more than a page is 

used to list the 34 particular numerical cases of the expression 


, 180° (T—2) 
ie T 


? 


when T has integer values from 8 to 36. 

Apart from the suggested change in the title, the above does not constitute 
a criticism of the book. The reviewer feels, however, that he cannot avoid mak- 
ing some criticism. He does this with the greatest reluctance because the book 
is a posthumous work, and therefore the author cannot answer any criticism 
with which he may have disagreed. 

After some interesting preliminary discussion, the first “mathematics” ap- 
pears on page 31. We quote here a whole paragraph. 

“An art structure in its formation (the process of being created) may be 
expressed as a series of sub-structures in their sequent development and ac- 
cumulation. 


=S,t, + (S, + Sy) t. + (Sy + Sy Syy)t, +°--- 


my QS 


ms (S, > Sy > Spy = + +++ + Sy) ty. 
S = substructure 
t, , t, , t, = consecutive moments in time.” 


Nowhere in the text is the meaning of the other symbols explained. On 
page 684, in-the glossary, compiled by Arnold Shaw, we read: “S used to denote 


*Italics by the author. 
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structural group in Positional Rotation* .... S also used to denote syn- 
chronization, or coordination of duration, or Symmetrization, the coordination 
of extensions, as in S(a:b).” This does not seem to be of much help! After some 
guessing, the reviewer is inclined to think that the left-hand side of the above 
expression means a summation from 1 to N, and should, in conventional terms, 


N 
read >’. The meaning of the symbol + is left for the reader to guess. If this 


a 
is a division sign, then what does it mean to divide one substructure by another? 
For that matter, what is the meaning of the product of a “structure” with time? 
The paragraph quoted above is followed by pages of purely verbal discussion. 
A mathematician is likely to find such a loose use of symbols rather irritating. 
As a further example we find the following on page 85. “A natural integer series is 
an infinite progression of all the integer numbers between one and any integer num- 


n 
ber approaching infinity. Y — 1, 2,3,---- 7.” Apparently in this case the symbol 


1 
» does not mean a summation; or does it? Again the glossary on page 684 fails 
to elucidate the point. Also on page 85 we read: “Continuity* is a finite portion 
. continuum” (sic!). Many more such peculiarities are found throughout the 
ook. 

The value and power of mathematics lies in its ability to provide general ex- 
pressions, which make it unnecessary to enumerate all particular cases. From 
this angle again we come to the conclusion that Schillinger’s book can hardly be 
called mathematical, although it contains page after page of numerical material. 

As we have already remarked, such a classificatory work may well prove to 
be useful in the future, not as a theory, which it decidedly is not, but as a sort 
of semi-empirical numerical material. Whether it will prove to be such or not, 
the future only will show. The amount of work done by the author in writing 
this book was tremendous. Possibly the author had some theoretical or practical 
applications in mind for the future. It is certainly to be regretted that his un- 
timely death leaves that conjecture unanswered forever. 


N. RASHEVSKY 
Committee on Mathematical Biology 
The University of Chicago 

*Italics by the author. 
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